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ABSTRACT 


The  final  phase  of  a  two  part  investigation  of  hypersonic  shock-wave 
boundary  layer  interaction  effects  from  a  secondary  shock  system  on  super¬ 
sonic/hypersonic  combustion  inlet  performance  has  been  completed.  In  support 
of  this  effort,  the  following  has  been  accomplishod : 

1.  The  computational  procedure  developed  in  Phase  I  for  estimating 
laminar  boundary  layer  properties  through  impinging  shock  induced 
laminar  boundary  layer  interactions  has  been  modified  to  account 
for  entropy  change  of  the  inviscid  streamlines  in  that  isentropic 
assumptions  were  used  in  the  initial  formulation.  Perturbation 
studies  on  the  effects  of  wall  shear  on  interaction  length  have 
identified  the  importance  of  the  profile  shape  factor  in  such 
calculations.  Initial  correlations  have  been  developed,  on  the 
basis  of  solutions  obtained  with  the  formulated  program,  for 
making  rapid  estimates  of  wall  pressure  distributions  through  the 
interaction  region  for  cases  where  the  laminar  boundary  layer 
remains  attached.  A  Users'  Manual  for  the  laminar  computational 
program  is  contained  herein. 

2.  Three  incident  shock  models  were  designed  for  use  in  the  experi¬ 
mental  program  conducted  in  the  Ames  3.5  Foot  Hypersonic  Tunnel. 
Tests  with  the  resulting  models  provided  laminar  and  turbulent 
boundary  layer  shock  Induced  separation  data  including  direct 
measurements  of  up  and  downstream  boundary  layer  profiles  as  well 
as  wall  static  pressure  distributions  through  the  interaction 
region.  Results  of  subsequent  data  analysis  were  compared  with 
applicable  theories  for  predicting  boundary  layer  separation 
characteristics. 

3.  Aided  by  analysis  of  the  eAperlmental  data,  a  simplified  flow  model 
was  formulated  for  the  case  of  a  two-dimensional  turbulent  boundary 
layer  incident  shock  Interaction.  Techniques  for  estimating  the 
change  in  boundary  layer  thickness  across  the  shock  impingement 
region  were  developed  which  give  results  In  general  agreement  with 
those  obtained  in  the  experimental  program.  Using  the  formulatod 
flow  model  co.  cept,  a  method  for  predicting  imminent  nopsrstlon 
was  developed  which  was  also  found  to  tie  consistent  with  experi¬ 
mental  results,  as  well  as  qualitatively  indicating  the  increased 
stability  of  the  turbulent  boundary  layer  at  decreased  wall 
temperature  and/or  Increased  upstream  power  law  profile  parameter 
values . 
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SECTION  I 


INTRODUCTION 

At  hypersonic  flight  speeds,  the  design  of  inlet  compression  surfaces 
necessarily  involves  the  use  of  secondary  shock  system  compression  to  achieve 
the  increased  combustor  entrance  pressure  and  temperature  levels  compatible 
with  high  levels  of  engine  internal  performance  at  low  values  of  external 
drag.  At  these  higher  flight  speeds,  viscous  effects  play  an  increasingly 
important  role  in  the  design  and  performance  of  inlet  component  compression 
surfaces.  While  various  methods  are  available  for  treating  the  resulting 
combined  viscous/inviscid  flow  fields  for  cases  of  continuous  adverse  pres¬ 
sure  gradient  resulting  from  single  family  wave  compression,  methods  for 
handling  viscous  effects  associated  with  oblique  shock  impingement  are  much 
more  limited  in  scope.  Treatment  of  such  problems  necessarily  involves 
consideration  of  boundary  layer  separation  phenomena.  At  supersonic  speeds, 
a  considerable  amount  of  adiabatic  wall  test  data  are  available  in  the  liters 
ture  which  can  be  used  in  assessing  the  probability  of  secondary  shock  in¬ 
duced  boundary  layer  separation  for  lower  speed  inlet  applications.  The 
applicability  of  such  lower  speed  techniques  to  higher  Mach  number  applica¬ 
tions,  however,  requires  evaluation  and  possibly  modification.  In  some 
instances  it  is  to  be  expected  that  new  methods  of  analysis  will  have  to  be 
formulated.  The  Investigations  described  in  this  report  have  been  addressed 
to  these  problt  i,  with  the  overall  objective  of  significantly  contributing 
to  their  resolution  as  well  as  providing  necessary  insight  to  further  re¬ 
quired  study.  In  accomplishing  this  purpose,  a  combined  analytical  and 
experimental  program  was  conducted  such  that  thecry/data  comparisons  could 
be  made  for  both  laminar  and  turbulent  boundary  layer-incident  shock  inter¬ 
actions  at  hypersonic  speeds.  The  computational  techniques  which  were  formu¬ 
lated  for  the  treatment  of  laminar  and  turbulent  interactions  were  found  to 
give  results  consistent  with  those  obtained  from  te3ts  conducted  with  the 
component  surface  models  at  Mach  7  and  10  as  well  as  with  data  reported  by 
other  investigators,  and  thus  are  considered  appropriate  for  current  engi¬ 
neering  use  in  the  design  and  performance  analysis  of  hypersonic  inlets. 
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SECTION  II 


ANALYTICAL  PROGRAM 

A.  LAMINAR  SHOCK/BOUNDARY  IAYER  ANALYSIS 

Reference  1  documented  the  formulation  of  a  method  for  the  analysir 
of  a  laminar  boundary  layer/incident  shock  wave  interaction  following  the 
initial  approach  of  Lees  and  Reeves  in  Reference  2.  A  User's  Manual  incor¬ 
porating  the  most  recent  modifications  to  the  resulting  program  is  presented 
in  Appendix  I  of  this  report.  The  calculation  procedure  simultaneously 
accounts  for  the  combined  Interaction  of  the  inviscid  and  viscous  flow  field 
and  establishes  the  extent  of  the  forward  propagation  of  the  upstream  pres¬ 
sure  disturbances  given  the  upstream  inviscid  conditions  and  the  shock  im¬ 
pingement  point.  Solutions  can  be  obtained  throughout  the  complete  laminar 
interaction.  The  results  of  such  calculations  have  generally  been  found  to 
agree  favorably  with  supersonic  test  data  as  indicated  in  Figure  1  of  thia 
report  for  the  data  of  Reference  3.  Similar  agreement  at  supersonic  test 
conditions  was  noted  earlier  in  Reference  1  with  the  data  of  Reference  4. 
Similar  comparisons  with  data  obtained  at  hypersonic  speeds  have  been  less 
successful;  due  in  part  to  the  onset  of  boundary  layer  transition  downstream 
of  the  interaction.  During  the  second  phase  of  this  study  attention  was 
directed  at  modifications  to  the  original  program  for  the  purpose  of  obtain¬ 
ing  better  theory/data  agreement  at  hypersonic  speeds  and  to  increase  the 
class  of  laminar  problems  for  which  complete  iterated  solutions  could  be 
obtained.  The  modified  computer  program  was  subsequently  used  for  hypersonic 
theory/data  comparisons  presented  in  Section  III-D  of  this  report. 

a.  Entropy  Rise  Effect 

Boundary  layer  growth  calculations  are  commonly  carried  out  on 
the  basis  of  lsentropic  external  flow  at  the  boundary  layer  edge.  Indeed, 
the  transformation  between  the  compressible  boundary  layer  equations  and  the 
low  speed  form  of  these  equations  upon  which  the  calculations  are  based  is 
dependent  upon  the  external  lsentropic  flow  assumption.  Although  this  re¬ 
striction  is  undesirable,  it  has  not  particularly  limited  the  engineering 
usefulness  of  the  results.  In  fact,  such  an  approximation  was  employed  in 
the  initial  formulation  of  the  Laminar  Shock  Boundary  Layer  interaction 
calculation  procedure  reported  in  Reference  1.  Even  at  hypersonic  speeds 
where  the  external  flow  may  be  highly  rotational  (due  to  a  curved  leading 
edge  shock  wave)  such  quantities  as  shear  stress  and  boundary  layer  momentum 
decrement  are  found  to  be  insensitive  to  the  external  total  pressure. 

Instead,  the  free  stream  stagnation  enthalpy  level  and  the  local  static 
pressure  are  found  to  be  the  principal  controlling  parameters. 

The  effect  of  total  pressure  change  in  a  shocn  boundary  layer 
Interaction,  however,  is  less  obvious  especially  at  hypersonic  speeds  and 
increased  shock  strengths.  At  low  Mach  numbers  the  laminar  boundary  layer 
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is  separated  by  a  rather  weak  incident  shock  and  thus  the  downstream  total 
pressure  is  usually  not  very  different  from  the  upstream  total  pressure; 
there  is  little  need  for  including  entropy  variation  effects  in  such  analy¬ 
ses.  At  higher  Mach  numbers,  however,  the  normal  shock  Mach  number  (of  the 
oblique  incident  shock  wave)  can  be  much  larger  and  therefore  the  total 
pressure  variation  is  of  concern;  particularly  at  increased  values  of  turn¬ 
ing  across  the  incident  shock  wave  system.  As  a  consequence  of  this,  a 
modification  has  been  made  to  the  Laminar -Shock/Boundary  Layer  Interaction 
computer  program  of  Reference  1  which  now  accounts  for  the  variation  in 
total  pressure  across  the  incident  shock. 

In  accomplishing  this  modification,  the  boundary  layer  edge  conditions 
are  still  based  upon  the  Prandtl-Meyer  equations: 

v  i  6  =  vQ  v  =  Prandtl-Meyer  angle  at 

boundary  layer  edge 

v  0  -  •'-  =  flow  angle  at  boundary 

layer  edge 

which  are  valid  upstream  and  downstream  of  the  incident  shock  wave  respec¬ 
tively.  The  value  of  vQ  is  determined  from  the  input  value  of  Mach  number; 
the  value  of  is  calculated  from  the  specified  shock  strength.  Specifi¬ 
cally,  the  boundary  layer  growth  and  edge  Mach  number  are  calculated  up  to 
2  in  Figure  2,  and  the  downstream  conditions  at  3  are  determined  by  a 
constant  pressure  boundary  shock  reflection  calculation.  A  knowledge  of  the 
change  in  v  and  6  across  the  shock  expansion  system  leads  to  the  following 
equation  for 


v"  =  V0  -  aes  -  <V2-V2A) 


where  A9g  is  the  given  flow  deflection  angle  through  the  shock  wave. 


The  calculation  procedure  has  been  formulated  so  that  values  of 
displacement  thickness  and  momentum  decrement  (momentum  decrement  equals 
2  * 

o  u  9  -  ,  where  0  is  a  boundary  layer  momentum  thickness)  are  held 

e  e 


invariant  together  with  static  pressure  from  2  to  3  in  Figure  2.  Other 
quantities  must  now  necessarily  be  discontinuous.  Whereas  the  isentropic 
calculation  procedure  of  Reference  1  permitted  continuity  across  the  inter¬ 
face,  the  entire  profile  within  the  boundary  layer  must  now  be  discontinuous 
as  well  as  the  boundary  layer  edge  values  of  total  pressure,  static  tempera¬ 
ture,  density,  Mach  number  and  velocity. 


Computed  results  using  both  the  original  isentropic  method  of  Reference  1 
and  this  subsequent  non-isentropic  modification  are  shown  in  Figure  3.  At 
these  conditions,  the  effect  on  the  static  pressure  distribution  of  account¬ 
ing  for  the  entropy  rise  agpears  to  be  small.  The  comparison  shown  is  for  a 
given  flow  deflection  of  5  across  the  impinging  shock  wave.  By  including 
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the  total  pressure  loss  (P^/P^  =  .923)  in  the  calculation,  tie  overall 
static  pressure  rise  has  been  reduced  somewhat  as  would  be  expected.  Although 
the  entropy  effect  would  be  more  noticeable  for  higher  impinging  shock 
strengths,  an  input  value  of  A0g  =  7.96  at  these  free  stream  conditions  ras 
found  to  he  too  strong  to  permit  a  complete  solution  to  be  obtained  with  the 
present  laminar  boundary  layer  model. 

b.  Strong  Shock  Convergence 

Modifications  to  the  Laminar-Shock/Boundary  Layer  Interaction  program 
have  been  made  to  increase  reliability  of  the  original  convergence  procedure. 
During  a  typical  calculation,  an  integration  of  the  boundary  layer/Prandt 1  - 
Meyer  equations  is  performed  from  a  downstream  point  (station  4)  to  the 
oblique  shock  incidence  point  (station  3).  This  integration  involves  assuming 
trial  values  of  momentum  thickness  6^  at  station  4  and  then  iterating  to 
obtain 


where  subscripts  2  and  3  denote  interface  values  immediately  upstream  and 
downstream  of  the  shock  wave  incidence  point.  Occasionally,  during  the 
solution  of  some  problems,  values  of  c*  may  be  selected  by  the  iteration 
procedure  which  cause  the  integration  to  fail  before  reaching  station  3. 

The  difficulty  is  associated  with  the  momentum  thickness  decreasing  to  zero 
as  shown  by  the  typical  example  of  Figure  4. 

When  this  situation  is  detected  in  the  program,  corrective  action  is  now 
taken  as  follows:  in  order  to  continue  the  iteration  for  the  appropriate 
64  a  negative  (fictitious)  value  of  momentum  thickness  is  calculated  by  the 
simple  relation 


Rcj/x  =  input  Reynolds  number  per  unit  length 

(See  Figure  4  for  definitions  of  X,  ,  and  X_„) 

A  plot  of  an  actual  iteration  process  using  this  procedure  is  presented  in 
Figure  5  exhibiting  the  present  capability  to  obtain  solutions  to  a  class 
of  problems  not  previously  possible.  It  should  be  pointed  out,  however, 
that  even  with  this  increased  capability,  complete  laminar  solutions  are 
not  always  possible  at  the  higher  shock  strengths.  This  analytical  result 
is  in  accord  with  the  experimental  results  which  conclusively  indicated  that 
increased  shock  strength  induced  transition  in  that  the  separating  laminar 
boundary  layer  subsequently  reattached  and  transitioned  to  turbulent  flow 
downstream  of  the  interaction. 
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c.  Perturbations  Studies 


To  better  understand  the  behavior  of  the  laminar  boundary  layer  under¬ 
going  a  sudden  pressure  rise,  it  is,  perhaps,  desirable  to  investigate  the 
computed  effects  resulting  from  the  perturbation  of  basic  parameters.  It 
was  felt  that  a  perturbation  study  would  demonstrate  the  significance  of 
certain  parameters  and,  in  particular  might  also  allow  determination  of  a 
possible  parametric  variation  which  would  result  in  a  smaller  interaction 
zone  with  hypersonic  flow.  The  results  of  this  study  are  discu.-'sed  in  the 
following  paragraphs. 

Primarily,  the  variation  of  the  shear  parameter  ^  has  been  studied. 
Modifications  were  thus  made  to  the  Laminar-Shock/Boundary  Layer  Interaction 
program  so  that  the  dimensionless  shear  parameter 


►A  U 

e 

could  be  arbitrarily  adjusted  by  the  relation 

^  -  Ax<.’  *  Ap  (4) 

where  V  is  the  unadjusted  "similar  solution”  value  of  t,  and  Ax  and  Ap  are 
arbitrary  constants  which  apply  over  the  region  downstream  (upstream  or  both 
upstream  and  downstream)  of  the  shock  wave  impingement.  Provision  was  also 
made  for  perturbation  of  the  kinetic  energy  dissipation  function 

*  =  Bx  '  BP  (5) 

With  such  an  approach,  however,  it  should  be  noted  that  a  condition 
imposed  on  all  of  the  shock  boundary  layer  calculations,  including  the 
current  perturbation  study,  is  that  at  the  beginning  and  at  the  end  of 
the  interaction  region  the  boundary  layer  must  reach  an  "equilibrium"  state. 
That  is,  the  dimensionless  velocity  profile  does  not  change  for  a  small 
distance  upstream  or  downstream  of  the  terminal  station.  This  implies  that 
the  KE/momentum  equation  (Equation  32  of  Reference  1)  must,  therefore  take 
on  the  value  of  zero  at  the  upstream  and  downstream  boundaries  since  dc/dX 
must  be  zero: 


2 

p  u  o 
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dG 

dX 


=  29  -  Gl  -  n 


(6) 


Since  it  is  also  usually  assumed  that  the  pressure  gradient  is  zero,  (n  =  0) 
at  the  terminal  stations,  the  above  equation  reduces  to 

20  -  Gl  z  0  (7) 
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Unfortunately  this  equation  will  not  remain  satisfied  if  one  of  the  parameters 
0,  G  or  t  is  varied  independently.  Furthermore,  at  the  upstream  and  downstream 
boundaries  of  the  interaction  region,  the  value  of  the  low  speed  form  factor, 

H  ,  is  chosen  so  that  Equation  (7)  is  satisfied.  This  requires  a  shift  in 
Htr  from  the  nominal  (flat  plate)  value  if,  for  example,  l  is  perturbed.  In 

the  first  three  sets  of  calculations  to  be  discussed,  the  value  of  t  was 
independently  perturbed  upward.  This  resulted  in  Htr  (and  6*)  taking  on  a 
larger  value.  In  the  fourth  set  of  calculations  both  L  and  0  were  selected 
so  that  other  profile  parameters  would  remain  constant.  For  the  four  sets 
of  perturbation  computations  carried  out,  the  results  are  shown  in  Figures  6 
through  9  and  are  discussed  below. 

1)  l  1.25t'  downstream  of  shock  impingement  only: 

5*  (i.e.  Htr)  increased  at  the  downstream  boundary  to  accommodate 
the  "equilibrium"  condition.  Consequently,  the  (initially 
unexpected)  result  was  that  the  shear  parameter  was  reduced  from 
the  intended  perturbed  value  to  a  level  which  duplicated  the 
unperturbed  value.  Just  downstream  of  the  shock  impingement,  the 
magnitude  of  the  negative  shear  was  found  to  increase  and  the 
pressure  gradient  was  slightly  steeper. 

2)  l  -  l'  .04  downstream  of  the  shock  impingement  only: 

This  resulted  in  a  somewhat  larger  perturbation  at  the  downstream 
boundary  from  the  previous  case.  Due  to  the  "equilibrium"'  condition, 
the  resulting  value  of  shear  actually  decreased.  6*  was  larger  over 
the  entire  region,  and  to  match  the  6*  at  the  shock  impingement 
point,  the  leading  edge  of  the  interaction  moved  forward. 

3)  t  ~  t  \  .04  both  upstream  and  downstream  of  the  shock  impingement: 
Similar  but  more  pronounced  effects  were  found  than  in  (2).  The 
length  of  the  interaction  was  increased  considerably. 

4)  l  =  V  .04,  .0315  downstream  of  the  shock  impingement: 

Using  Equation  (7)  perturbation  values  of  t  and  0  were  selected  so 
that  the  other  profile  parameters  (G,  Htr,  n,  etc.)  would  remain 
unchanged  (at  the  downstream  ooundary)  .  In  this  case  a  truly  in¬ 
creased  value  of  shear  was  obtained  over  the  entire  range  of  the 
perturbation.  However,  the  pressure  distribution  and  the  length 
of  the  interaction  were  only  slightly  affected. 

It  is  concluded  from  this  study  that  the  parameters  l  and  0  (and 
probably  G  and  e/3  also)  play  second  order  effects  compared  to  the  form  factor, 
Htr.  Interaction  effects  with  the  free  stream  are  determined  by  the  rate  of 
boundary  layer  growth.  As  the  value  of  momentum  thickness  does  not  change 
drastically  over  the  interaction  region,  changes  in  the  value  of  6*  follow 
closely  changes  in  the  parameter  Htr.  For  the  turbulent  boundary  layer, 
the  form  factor  is  smaller  than  for  the  laminar  layer.  The  corresponding 
small  interaction  region  of  the  turbulent  boundary  layer  as  observed  in  a 
subsequent  section  is  consistent  with  these  laminar  perturbation  solution 
results  in  that  the  length  of  the  interaction  seems  to  be  primarily 
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dependent  upon  the  change  in  layer  thickness  or  more  specifically,  on  H 

t  r 

d.  Interaction  Correlation 

Bv  making  certain  approximations,  in  Reference  7  Erdos  and  Pallone 
developed  a  correlation  of  the  axial  pressure  distribution  for  a  "free 
interaction"  boundary  layer  separation  region.  Such  a  "free  interaction"  will 
exist  upstream  of  an  impinging  shock  wave  if  its  strength  is  sufficient  to 
separate  the  boundary  layer  and  form  the  familiar  pressure  plateau.  Unfortu¬ 
nately,  the  extent  of  forward  feed  of  the  pressure  rise  (or  "free  interaction") 
has  not  been  correlated  successfully.  This  forward  feed,  or  conversely,  the 
degree  to  which  the  pressure  plateau  is  evident  is  a  function  of  the  shock 
wave  strength.  By  introducing  a  shock  strength  correlation  parameter,  the 
concepts  of  Erdos  and  Pallone  are  here  extended  to  include  both  the  separation 
and  reattachment  zones  in  the  shock-boundary  layer  interaction.  Thus  the 
axial  pressure  distribution,  including  the  extent  of  upstream  and  downstream 
separation,  is  correlated  for  various  values  of  the  shock  strength  parameter. 
Based  on  laminar  boundary  layer  calculations,  the  correlation  appears  to  be 
valid  for  a  limited  range  of  conditions. 

The  approximations  used  in  the  correlation  development  are  presented 
below : 


1)  The  normal  gradient  of  shear  stress  at  the  wall,  throughout  the 
interaction  region,  is  related  to  the  wall  shear  and  boundary 
layer  thickness  at  some  reference  point,  and  a  universal  function  f2« 


(8) 

The  subscript  'r'  denotes  the  reference  point,  which  Erdos  and 
Pallone  take  as  the  beginning  of  the  pressure  rise.  Most  signifi¬ 
cant,  the  function  f„  is  assumed  to  apply  at  all  Reynolds  number- 
Mach  number  conditions;  AX  is  a  characteristic  length  of  the  inter¬ 
action  region. 


2)  Similarly,  the  boundary  layer  edge  slope  is  related  to  a  reference 
value  of  5*  and  a  distance  from  the  reference  point  by  a  second 
universal  function,  f^ 


(9) 
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is  approximated  by 


v 
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Approximations  (9)  and  (10)  imply  that  the  length  of  the  interaction  is  small 
(or  constant)  relative  to  distance  from  the  effective  origin  of  the  boundary 
layer.  That  is,  the  change  in  boundary  layer  edge  slope  and  t  /p  due  to  the 
natural  boundary  layer  growth  (without  a  shock  wave  impingement!)  must  be 
small  compared  to  the  changes  imposed  by  the  sudden  pressure  rise.  The 
turbulent  boundary  layer  meets  these  requirements  more  closely  than  the 
laminar  boundary  layer  since  d^*/dX  and  Tw/p  nominally  vary  with  the  .2 
power  of  X  rather  than  the  .5  power.  Approximation  (11)  limits  the  correla¬ 
tion  to  weak  interactions. 


In  spite  of  possible  limitations  of  the  correlation  formula,  especially 
for  a  laminar  interaction,  it  is  worthwhile  to  explore  the  calculated  results 
using  the  Laminar-Shock/Boundary  Layer  Interaction  theory.  In  order  to 
generalize  Erdos  and  Pallone's  work  to  the  entire  interaction  region,  a 
shock  strength  parameter 


(12) 


ABg  -  flow  angle  deflection 
through  the  impinging  shock 
wave,  in  radians 


is  introduced.  The  pressure  and  distance  parameters  are: 


q 
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P  u 
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(13) 


(14) 


Erdos  and  Pallone  chose  the  reference  point,  denoted  by  'R',  as  the  beginning 
of  the  interaction  region.  However,  because  of  the  difficulty  of  locating 
the  beginning  of  the  smooth  pressure  rise,  it  was  elected  to  calculate  the 
reference  values,  and  ,  at  the  shock  impingement  point  assuming 

undisturbed  flow,  instead.  R 


The  resulting  correlation  is  found  to  be  only  moderately  satisfactory 
for  laminar  layerr.  Calculated  results  for  a  shock  strength  parameter  of 
unity  and  a  range  of  Mach  number  and  Reynolds  number  conditions  are  shown 
in  Figures  10  and  11.  Figure  10  shows  the  actual  calculated  wall  pressure 
distributions  for  an  oblique  shock  wave  impinging  on  a  flat  plate.  Note 
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that  the  impinging  shock  is  quite  weak  so  that  in  these  examples  no  pressure 
plateau  is  formed.  In  Figure  11  the  same  data  are  plotted  in  terms  of  the 
correlation  parameters.  The  MQ  z  10,  Reyg^  =  1.6  X  10^  condition  is  in¬ 
cluded  for  comparison  although  Approximation  (11)  is  badly  violated,  in  this 
particular  case. 

It  appears  that  the  correlations  obtained  have  applicability  in 
regions  of  low  Mach  number  and  high  Reynolds  number  for  a  laminar  boundary 
layer  and  thus  in  some  instances  could  be  used  in  place  of  more  formal 
computer  solutions. 

n.  turbulent-shock/bounda RY  LAYER  SEPARATION  FORMUUTION 

Comparison  of  schlieren  and  wall  static  pressure  data  (Section  III  D  a, 
Figures  75  through  80)  document  contrasting  differences  between  the  laminar 
and  turbulent-incident  shock  wave/boundary  layer  interactions.  Of  particular 
importance  in  the  formulation  of  simplified  flow  models  is  tne  observation 
that  the  shock  impingement  and  reflection  from  a  non-separating  laminar 
boundary  layer  occurs  more  or  less  as  a  constant  pressure  boundary  phenomena, 
whereas  with  the  turbulent  interaction  the  impinging  shock  penetrates  deep 
into  the  viscous  flow  and  reflects  in  a  manner  similar  to  that  associated 
with  a  non-uniform  flow  field  interaction.  Such  differences,  coupled  with 
the  fact  that  data  obtained  under  this  program,  as  reported  in  Section  III-D, 
have  consistently  documented  the  stability  of  the  non-ad iaba t ic  wall  hyper¬ 
sonic  turbulent  boundary  layer  over  a  range  of  shock  impingement  strengths 
such  as  might  be  associated  with  hypersonic  inlet  operation,  led  to  the 
formulation  of  the  simplified  flow  model  depicted  in  the  schematic  of 
Figure  12. 

With  respect  to  analysis  in  connection  with  the  simplified  two- 
dimensional  flow  model  of  Figure  12,  the  following  key  assumptions  are  made: 

1.  The  inviscid  and  viscid  conditions  up  to  a  are  known. 

2.  There  is  no  significant  mass  addition  in  the  region  a  to  B. 

3.  Wall  friction  is  negligible  in  the  region  A  to  B  in  regard  to 
its  influence  on  profile  shape. 

4 .  Boundary  layer  edge  conditions  and  flow  direction  can  be  deter 
mined  from  the  known  incident  shock  strength. 

5.  The  reflected  shock  strength  is  that  required  to  turn  the  boundary 
laver  edge  flow  back  parallel  to  at  B. 

6.  The  wall  temperature  to  stream  total  temperature  is  constant. 

7.  Constant  pressure  mass  addition  occurs  in  the  region  B  to  C. 
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8.  The  distorted  boundary  layer  profile  at  B  mixes  out  to  a  typical 
turbulent  profile  at  C. 

In  spite  of  the  fact  that  such  a  formulation  is  an  over-simplification 
of  this  complex  problem,  it  is  believed  that  results  obtained  from  subsequent 
related  analysis  will  be  of  engineering  usefulness  in  current  hypersonic 
inlet  design  applications  until  more  detailed  analyses  are  conducted. 

Many  of  the  preceding  flow  model  assumptions  are  an  outgrowth  of  data 
analysis  subsequently  presented  in  Section  III-D  of  this  report.  For 
instance,  in  Figure  13  it  is  seen  that  such  a  flow  model  leads  to  reasonable 
approximations  for  the  inviscid  boundary  interaction  length  while  providing 
indirect  substantiation  for  assumptions  4  and  5.  The  degree  to  which 
assumptions  7  and  8  are  experimentally  justified,  beyond  expediency  of 
analysis,  can  be  assessed  from  the  wall  and  boundary  layer  edge  pressure 
distribution  comparisons  of  Figure  14  and  probe  station  3  boundary  layer 
profile  measurements  reported  and  discussed  in  Section  III-D.  The  schlieren 
photographs  of  Figures  15  through  17  document  visual  similarities  with  the 
simplified  flow  model  of  Figure  12  and  suggest  that  such  a  formulation  might 
be  acceptable  for  incorporation  with  combined  Method  of  Character ist ics- 
Boundary  Layer  analysis  of  the  inlet  flow  field,  providing  boundary  layer 
integral  properties  (particularly  displacement  thickness)  can  be  reasonably 
estimated  immediately  downstream  of  the  interaction. 

Consequently,  during  the  course  of  this  study,  two  methods  of  analysis 
constrained  by  the  preceding  assumptions  were  developed  for  estimating 
boundary  layer  integral  property  change  across  the  immediate  region  of  an 
impinging  shock  interaction.  Numerical  results  obtained  from  both  methods 
are  similar  and  agree  favorably  with  data  obtained  under  this  program. 

Limited  checks  with  data  reported  by  other  investigators  indicated  similarly 
favorable  agreement. 

In  the  first  method,  it  is  assumed  that  compresrtble  flow  profiles  can 
be  adequately  represented  by  the  modified  Crocco  enthalpy  distribution  and 
the  power  law  velocity  expression: 

u 
u 
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Momentum  and  continuity  are  balanced  across  the  boundary  layer  control 
volume  from  a  to  B  in  Figure  12.  This  is  an  iterative  calculation  to 


determine  the  downstream  shape 
solution  of :  ^ 

CONTINUITY :  *3 


factor  at  B  in  Figure  12  through  simultaneous 
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The  second  method,  described  in  detail  in  Reference  5  is  cn  integral 
method  approach.  The  incompressible  moment  of  momentum  equation  is  solved 
to  obtain  an  expression  for  the  change  in  incompressible  shape  factor. 

Simple  power  law  profiles  in  the  incompressible  plane  are  transformed  using 
the  Stewartson  transformation  to  obtain  a  compressible  solution.  The  Crocco 
temperature  distribution  is  utilized  to  obtain  non-adiabatic  wall  solutions. 
The  latter  indicate  a  significant  increase  in  stability  (in  terms  of  shock 
strength)  due  to  wall  cooling.  The  computational  procedure  as  prepared  for 
desk  calculation  is  presented  in  Appendix  II  along  with  the  necessary  integral 
computational  curves. 

Numerical  results  are  compared  with  incident  shock  Model  C  data  in 
Figure  18.  The  shock  pressure  ratio  and  Mach  number  ratios  used  in  the 
analytical  estimates  are  based  on  oblique  shock  solutions  for  the  shock 
generating  cowl  angle  of  attack.  The  upstream  boundary  layer  profile 
characteristics  used  in  the  analysis  are  based  on  results  obtained  from 
single  probe  profile  traverses.  The  experimentally  determined  boundary 
layer  thickness  ratios  were  established  through  interpretation  of  schlieren 
photographs  supported  by  the  upstream  profile  measurements.  Considering  the 
simplifying  assumptions  inherent  in  the  analytical  model  upon  which  the 
formulated  calculation  procedure  is  based,  the  degree  of  agreement  indicated 
between  the  numerical  results  and  the  data  over  such  a  wide  range  of  shock 
strengths  is  considered  most  gratifying. 

The  trend  of  the  data  and  the  numerical  results  exhibited  in  Figure  18 
suggest  an  interesting  trade-off  occurs  between  boundary  layer  thinning  due 
to  compression  of  the  boundary  layer  "stream  tube"  and  thickening  due  to 
profile  distortion.  At  the  weaker  shock  strengths,  there  is  little  profile 
distortion  and  the  downstream  boundary  layer  thins  as  a  result  of  the 
continuity  balance  associated  with  the  compression  process.  As  the  shock 
strength  is  increased,  the  downstream  profile  becomes  significantly  enough 
distorted  to  promote  a  thickening  trend.  Ultimately,  as  separated  flow 
conditions  are  approached  the  profile  distortion  dominates,  with  the  net 
result  that  the  boundary  layer  thickness  ratio  begins  to  increase  as  the 
shock  strength  is  increased. 

If  it  is  tacitly  assumed  that  boundary  layer  separation  is  imminent 
upon  distortion  of  the  incompressible  profile  at  B  (Figure  12)  to  a  value 
of  N  '  1,  then  a  limiting  Mach  number  ratio  for  attached  boundary  layer 
across  the  shock  interaction  can  be  computed  using  the  procedures  outlined 
in  Appendix  II  with  results  as  presented  in  Figure  19.  (The  choice  of  a 
limiting  value  of  N  =  1  results  from  the  observation  that  reverse  profile 
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curvature,  similar  to  that  encountered  as  separated  flow  conditions  are 
approached,  occurs  for  N  <  1).  In  examining  the  curves  of  Figure  19,  it  is 
noted  that  the  predicted  influence  of  wall  temperature  ratio  is  in  accord 
with  experimental  observations.  At  supersonic  speeds,  where  near  adiabatic 
wall  temperatures  are  generally  the  case,  less  compression  (in  terms  of  Mach 
number  reduction)  is  predicted  before  boundary  layer  separation  is  imminent. 
On  the  other  hand,  at  hypersonic  speeds,  where  wall  cooling  is  a  physical 
necessity,  and  particularly  at  the  low  wall  temperature  ratios  associated 
with  shock  tunnel  testing,  significantly  more  compression  potential  is 
predicted  before  the  on-set  of  boundary  layer  separation. 


Accepting  the  results  of  Figure  19,  it  is  then  desirable  to  re-cast  such 
separation  criteria  into  a  form  perhaps  more  suitable  for  general  application. 
Boundary  layer  separation  pressure  plateau  levels  are  frequently  expressed 
as  a  function  of  parameters  such  as  Mach  number,  skin  friction  coefficient, 
and  Reynolds  number  just  upstream  of  the  shock  impingement  point.  (See 
for  instance,  Reference  6  and  7.)  On  the  basis  of  the  experimental  results 
reported  in  Section  III  D  however,  it  is  apparent  thit  a  wide  range  of 
non-separa ted  boundary  layer  shock  induced  pressure  rise  ratios,  such  as 
might  be  employed  in  hypersonic  inlet  compression  processes,  are  attainable 
before  a  turbulent  separation  pressure  plateau  occurs.  Consequently,  it 
seems  appropriate,  and  is  perhaps  of  greater  engineering  usefulness  to 
express  such  separation  criteria  in  terms  of  the  inviscid  shock  strength 
for  which  separation  is  imminent  at  a  prescribed  level  of  upstream  wall 
shear.  In  accomplishing  this  end,  Figures  2U  through  22  were  prepared  by 
combining  the  separation  prediction  of  Figure  19  with  the  turbulent  flat  plate 
skin  friction  law  developed  by  Spence  in  Reference  8  which  can  be  expressed 
for  compressible  flow  as: 
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(See  page  32  of  Section  III-D-a  lor  additional  information.)  This  particular 
skin  friction  law  was  selected  because  of  its  generality  with  respect  to  N 
and  because  of  the  variation  of  N  encountered  in  the  test  results  of  this 
program.  The  predicted  sensitivity  of  turbulent  boundary  layer  separation 
to  N  (Figure  19)  amplifies,  for  the  present  application,  the  relative 
scatter  in  correlation  plots  of  this  power  law  parameter  as  a  function  of 
Reynolds  number  (curves  presented  in  References  9  and  10  for  example)  .  It 
was  thus  decided  best,  at  present,  not  to  reflect  a  correspondence  between 
the  power  law  profile  index  and  Reynolds  number  in  separation  curves  such 
as  those  of  Figures  20  to  22. 


In  using  the  latter  curves  to  assess  the  probability  of  incident  shock 
induced  turbulent  boundary  layer  separation,  it  should  be  noted  that  for  any 
given  impinging  shock  strength  (/9)  and  upstream  flow  conditions,  the 
corresponding  values  of  2  given  in  Figures  20  to  22  represent,  on 

C  (Reb)  N+l 
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the  basis  of  the  formulated  analysis, 
separation.  Thus  for  conditions  with 


the  minimum  values  required  to  avoid 
greater  2  no  separation 
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is  predicted,  and  correspondingly,  at  lower  values  of  C^.(Re  )  N+l  j  separation 
is  presumed  imminent.  This  criteria  used  in  Section  III-D-a  in  '■onjunction 
with  shock  impingement  data  obtain'd  under  this  program,  and  is  also  compared 
against  data  obtained  by  other  investigators  (Figure  101)  wherein  required 
upstream  boundary  laver  profile  ata  as  well  as  shock  strengths  required  for 
separation  were  reported.  The  rel  five  agreement  between  the  experimental 
results  and  the  predictions  are  considered  to  substantiate  the  engineering 
usefulness  of  this  relatively  simple  formulation. 


SECTION  III 
EXPERIMENTAL  PROGRAM 


A.  MODE  IS  AND  INSTRUMENTATION 


a .  Mode  1  C 

(1)  Design  and  Fabrication 

Model  C  is  basically  a  twice  scale  version  of  Model  A  which  was  shock 
tunnel  tested  earlier  in  the  program  at  Cornell  Aeronautical  Laboratory  (CAL) 
and  at  General  Electric  Missiles  and  Space  Division  (MSD)  as  reported  in 
Reference  1.  For  the  Model  A  shock  tunnel  tests  emphasis  was  placed  on 
static  pressure,  heat  transfer,  and  skin  friction  distributions  for  shock 
induced  boundary  layer  interaction  phenomena.  The  emphasis  during  Model  C 
testing  in  the  NASA-Ames  3.5  foot  Hypersonic  Wind  Tunnel  was  on  boundary 
layer  profile  and  ramp  pressure  measurements. 

Model  C  was  designed  and  fabricated  as  two  independent  major  assemblies: 
the  fixed  3  degree  ramp  assembly,  Figures  23  and  2d,  and  the  remotely  position- 
able  flat  plate  cowl  assembly,  Figures  25  and  26. 

The  fixed  3  degree  ramp  assembly  consists  of  an  18"  wide  X  48"  long 
7075  Aluminum  ramp  plate  (forward  portion  ^atercooled)  with  a  detachable 
tool  steel,  water  cooled,  sharp  leading  edge  component,  and  a  detachable 
cast  aluminum  undercarriage  housing  with  a  steel  U  section  adapter.  Top 
and  bottom  view  photographs  of  the  ramp  assembly  (before  instrumentation) 
are  shown  in  Figures  23  and  24.  The  aluminum  undercarriage  housing  rigidly 
supported  the  ramp  plate  and  housed  the  probe  actuators  and  other  attendant 
instrumentation.  Final  finishing  of  the  ramp  aerodynamic  surface  was 
accomplished  with  the  undercarriage  housing  installed.  The  transition  from 
the  steel  leading  edge  component  surface  to  the  aluminum  ramp  surface  was 
perfectly  smooth  and  remained  so  during  the  79  test  runs.  The  U  shaped 
steel  adapter,  attached  to  the  undercarriage  housing  and  the  bottom  of  tie 
ramp  plate,  mounted  on  the  steel  NASA-Ames  Pressure  Transducer  Pod  Housing 
which  in  turn  was  supported  by  the  tunnel  strut  system 
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The  cowi  assembly  consists  of  an  18"  wide  X  4l"  long  7075  Aluminum 
cowl  plate  (forward  portion  watercooled)  with  a  detachable  tool  steel, 
watercooled,  sharp  leading  edge  component,  a  7  long  aluminum  cowl  extension 
plate,  cowl  side  plates,  a  4130  Steel  support  and  pivot  assembly,  and  a 
cowl  drive  assembly.  Fop  and  bottom  view  photographs  of  the  cowl  plate 
with  cowl  extension  plate  and  support  fittings  installed  are  shown  in  Figures 
25  and  26. 

(2)  Instrumentation 

The  Model  C  ramp  was  instrumented  with  69  static  pressure  taps  and 
11  surface-embedded  thermocouples.  The  location  and  identification  of  the 
taps  and  thermocouples  are  shown  in  Figure  27.  The  .063  OD  X  .042  ID 
leadout  tubing  installation  and  one  probe  actuator  installation  are  shown 
in  the  photograph  of  Figure  28.  All  the  .063  0D  pressure  tubes  were  press- 
fitted  into  pre-drilled  holes  in  the  ramp  plate  and  honed  flush  with  the 
top  surface.  Each  tube  was  painted  with  Glvptol  where  it  emerged  from  the 
underside  of  the  ramp  plate  to  insure  an  airtight  seal. 

The  Model  C  cowl  rlate  was  instrumented  with  10  static  pressure  taps 
and  4  surface-embedded  thermocouples.  The  location  and  identification  of 
the  taps  and  thermocouples  are  shown  in  Figure  29. 

The  small  probe  actuator  used  to  drive  the  NASA-Ames  pitot  pressure- 
static  pressure-total  temperature  combination  probe  is  shown  in  the  photograph 
of  Figure  30  (shown  approximately  full  size).  The  actuator  consists  of  a 
17,100  ItPM  DC  motor  and  appropriate  gear  reduction  to  provide  .80  inch 
travel  in  approximately  10  seconds.  This  rapid  travel  time  is  of  the  same 
order  as  the  probe  instrumentation  response  time.  Thus  the  boundary  layer 
pitot  pressure  survey  data  points  were  obtained  shortly  after  the  probe 
instrumentation  equilibrated.  The  actuator  included  a  10  turn  potentiometer 
for  measuring  probe  position.  The  potentiometer  was  connected  as  two  arms 
of  a  bridge  circuit  which  enabled  zeroing  at  any  position.  As  the  probe 
height  changed,  the  unbalanced  bridge  voltage  was  displayed  on  a  digital 
voltmeter  to  provide  continuous  probe  position  readings.  Pre-test 
calibration  of  the  probe  height  positioning  circuits  indicated  that  the 
probe  could  be  repositioned  at  any  height  with  an  accuracy  of  ^-.002". 

Also  provided  on  the  probe  actuator  were  two  limit  switches.  These  switches 
limited  the  probe  travel  so  that  the  probe  would  not  emerge  from  the  surface 
of  the  ramp  beyond  the  design  limit  of  .7"  and  would  not  become  jammed 
against  the  ramp  surface  during  probe  withdrawal. 

A  pitot  pressure-static  pressure  total  temperature  combination  probe 
developed  by  NASA -Ames  was  used  for  the  ramp  boundary  layer  profile  surveys. 
The  three-prong  probe  is  shown  in  Figure  31.  Views  of  the  probe  and  a 
typical  installation  on  the  ramp  surface  are  depicted  in  the  photographs  of 
Figure  32.  Alternate  static  probe  y"  was  used  for  all  data  presented  in 
this  report.  The  combination  probe  was  remotely  positioned  at  various 
points  across  the  boundary  layer  and  its  height  position  and  the  three 
measured  quantities  were  simultaneously  recorded  on  the  Beckman  recording 
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system.  The  output  tape  from  the  Beckman  system  was  then  processed  using  the 
NASA  data  acquisition  program.  The  probe's  height  position,  pitot  pressure, 
static  pressure  and  total  temperature  wc-i  e  printed  out  for  use  in  subsequent 
data  reduction. 

(3)  Boundary  Layer  Trips 

For  the  Model  C  -  NASA  Ames  tests  several  factors  known  to  influence 
transition  were  fixed  by  the  3.5  foot  tunnel  test  facility  set-up.  These 
were  Mach  number,  wall  to  free  stream  temperature  ratio,  and  free  stream 
turbulence  le-'el.  The  magnitude  of  the  Reynolds  number,  based  on  distance 
from  the  model  leading  edge,  was  limited  by  total  conditions  available  and 
model  size. 

To  promote  boundary  layer  transition,  a  series  of  boundary  layer  trip 
configurations  were  utilized.  Limited  testing  was  accomplished  on  several 
trip  configurations  as  indicated  in  the  schematic  sketch  of  Figure  33. 

However,  trip  configuration  #5  (which  incorporates  trip  configuration  #3) 
proved  most  effective  and  was  used  for  the  balance  of  the  Mach  10  testing 
program . 

Trip  configuration  #3  consisted  of  l/8"  diameter  steel  balls  on  3/8" 
centers  being  embedded  into  and  then  high  temperature  brazed  onto  a  .020" 
thick  steel  strip.  The  overall  height  of  the  trip  was  .130".  This  configura¬ 
tion  was  secured  to  the  ramp  surface  with  epoxy  and  also  at  12  intermediate 
positions  across  the  width  of  the  ramp  with  small  screws. 

Trip  configuration  #5  was  comprised  of  Trip  configuration  #3  (at  ramp 
station  X  =  9.8")  plus  two  rows  (at  ramp  stations  X  =  4"  and  X  =  7")  of 
#2-64  Fillister  Head  screws  on  3/8 "  centers  across  the  width  of  the  ramp. 

(4)  Installation 

The  Model  C  ramp  assenbly,  without  the  cowl  assembly,  was  installed  in 
the  tunnel  for  the  first  34  test  runs.  The  installation  is  shown  in  the 
photograph  of  Figure  34.  The  ramp  was  axially  oriented  in  the  test  section 
so  as  to  utilize  three  of  the  five  available  schlieren  windows.  The  three 
windows,  designated  windows  #2,  #3,  and  #4,  are  centered  at  test  section 
stations  46",  63",  and  80"  respectively.  The  leading  edge  of  the  ramp  was 
located  at  test  section  station  41.25"  or  approximately  l"  downstream  of  the 
upstream  edge  of  the  #2  window  glass.  A  special  steel  panel  was  fitted  into 
the  #5  window  frame  (centered  at  station  102")  for  instrumentation,  accessory, 
and  watercooling  lines. 

Prior  to  the  ramp  installation,  the  NASA-Ames  Pressure  Transducer  Pod 
Housing  (with  provisions  for  72  Statham  pressure  transducers  with  electrical 
disconnect  assemblies)  was  mounted  on  the  ramp  undercarriage  adapter  and 
then  all  model  to  transducer  instrumentation  hookups  were  completed  and 
checked  out. 
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The  Model  C  cowl  assembly  was  installed  in  the  tunnel  for  test  runs 
35  through  66.  The  Model  C  ramp  and  cowl  installation  is  shown  in  the 
photograph  of  Figure  35.  Complete  Model  C  assembly  and  installation  details 
are  documented  in  the  reduced  size  assembly  drawing  of  Figure  36. 

A  special  4130  Steel  hatch  plate,  lo"  wide  X  40 "  long  X  2  l/4"  thick, 
was  fabricated  as  part  of  the  cowl  assembly.  This  plate  fitted  into  the 
existing  hatch  opening  on  top  of  the  tunnel  test  section  and  served  as  the 
mounting  plate  for  the  cowl  support  yoke  and  the  cowl  drive  assembly. 
Vertical  adjustment  (3/4")  of  the  cowl  was  provided  by  spacers  inserted 
between  the  hatch  plate  and  the  cowl  support  yoke.  Axial  adjustment  (l" 
upstream)  of  the  cowl  was  also  provided. 

The  cowl  was  driven  through  13  degrees  angle  of  attack  by  means  of  a 
commercially  available  worm  gear  jack  assembly  (2  ton  capacity,  6"  travel) 
which,  in  turn,  was  driven  by  a  direct-current,  variable  speed  motor.  Both 
the  jack  assembly  and  the  motor  assembly  were  mounted  on  the  hatch  plate 
outside  the  tunnel.  The  screw  of  the  jack  assembly  was  attached  to  a  1.5L)" 
diameter  4140  Steel  shaft  extending  through  the  bearing  surface  and  "0" 
ring  pressure  seal  in  the  hatch  plate.  The  shaft  in  turn  was  bolted  and 
pinned  to  a  steel  yoke  which  carried  a  cam  follower.  The  follower  rode 
in  the  cam-slot  of  an  inverted  T-soction  (slide  follower)  bolted  to  the 
forward  and  outboard  side  of  the  cowl  plate.  The  cowl  angle  position  was 
determined  by  a  linear  motion  potentiometer  mounted  outside  the  tunnel  on 
the  jack  assembly  support  housing.  The  shaft  of  the  potentiometer  was, 
in  turn,  connected  to  the  cowl  drive  shaft.  The  resolution  of  the 
potentiometer  provided  cowl  angle  position  to  within  0.1  degree. 

A  sheet  metal  streamlined  heat  shield  (l/8"  thick  4130  steel)  was 
attached  to  the  underside  surface  of  the  hatch  plate  to  protect  the  outboard 
portion  of  the  drive  assembly  strut  from  heat  distortion.  A  similar  shield 
was  attached  to  the  cam-follower  yoke  to  protect  the  inboard  portion  of  the 
strut  and  the  cam-follower  yoke  assembly.  The  two  shields  overlapped  with 
sufficient  clearance  for  the  required  4  1/2"  movement  of  the  strut.  The 
flexible  hose  water  lines  leading  from  the  cowl  were  routed  to  fittings 
threaded  into  the  underside  of  the  hatch  plate. 

No  major  problems  were  encountered  in  the  installation  and  operation 
of  the  cowl  assembly.  However,  one  minor  modification  was  required  after 
the  first  cowl-on  test  run.  A  slight  deflection  of  the  cowl  assembly  drive 
shaft  caused  the  cam  follower  to  move  in  the  cam  slot,  allowing  an  angular 
shafting  of  the  cowl.  This  problem  was  solved  by  restraining  the  shaft 
with  a  brass  bushing  supported  within  the  outboard  heat  shield. 

The  Model  C  ramp  and  cowl  assembly  remained  steady  during  the  79  test 
runs;  no  vibration  or  shifting  of  the  model  was  evident.  Both  the  3.5  foot 
tunnel  and  the  model  started  without  difficulty  in  spite  of  a  frontal  area 
blockage  ratio  of  17%. 
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b.  Model  D 


Model  D  was  a  modified  version  of  Model  C  and  consisted  of  two  60 
degree  swept  vane  compression  surfaces  mounted  on  the  fixed  3  degree  ramp 
surface  of  the  Model  C  ramp  assembly  as  shown  in  Figure  37.  The  purpose  of 
testing  Model  D  was  to  investigate  the  interaction  between  vane  generated 
shocks  and  the  ramp  boundary  layer  such  as  exists  in  three-dimensional 
SCRAMJET  inlet  configurations. 

The  vanes  are  designated  Right  Vane  and  Left  Vane  as  viewed  from  the 
aft  end  of  the  model  looking  upstream.  The  vanes  were  manually  pivoted  from 
6  degrees  to  8  detrees  on  the  ramp  surface  (and  locked  in  position)  to  provide 
variable  strength  vane  shock  pressure  rises.  Boundary  layer  trip  configura¬ 
tion  # 5  (Figure  33)  is  shown  installed  on  the  ranp  surface  (top  photograph) 
for  one  of  the  test  runs.  A  total  of  11  test  runs  (runs  24  through  34)  were 
made  with  the  Model  D  configuration  prior  to  the  installation  of  the  Model  C 
cowl  assembly. 

Vane  geometry  and  instrumentation  location  and  identification  are 
documented  in  the  schematic  sketch  of  Figure  38.  The  internal  wedge  angle 
of  each  vane  was  6  degrees  when  in  the  design  position  as  shown  in  the 
installation  configuration  of  the  schematic  sketch.  The  sharp-pointed 
leading  edge  of  each  vane  was  located  16.65  inches  from  the  ramp  leading 
edge  with  the  vane  set  in  the  design  p  -,ition.  The  right  vane  was  instrumented 
with  17  static  pressure  taps  on  the  cc  ^ression  surface  and  2  static  pressure 
taps  on  the  top  surface.  The  left  vane  was  instrumented  with  5  static  pressure 
taps  located  in  the  same  positions  as  5  of  the  taps  in  the  right  vane  and  was 
used  to  check  symmetry  of  vane  internal  flow  field.  Five  surface-embedded 
thermocouples  were  also  installed  in  the  left  vane  in  the  same  positions  as 
5  of  the  static  pressure  taps  on  the  right  vane.  These  thermocouples  were 
installed  to  monitor  the  vane  temperatures  so  as  to  prevent  overheating  of 
the  several  thin  sections  of  the  7075  Aluminum  vanes. 

The  outboard  sides  of  the  vanes  were  milled  out  as  shown  in  Figure  39 
to  facilitate  drilling  the  holes  for  the  pressure  and  temperature  instrumenta¬ 
tion.  The  pressure  tubing  was  press  fitted  into  the  pre-drilled  holes  in 
the  compression  surface  and  honed  flush.  The  right  vane  is  shown  fully 
instrumented  with  the  pressure  tubing  held  secure  with  nichrome  strips 
tack  welded  to  steel  strips  which  were  in  turn  bolted  to  the  aluminum  vane. 

c.  Model  E-3 

(1)  Design  and  Fabrication 

Model  E-3  is  a  two  dimensional  isentropic  compression  ramp  with  a  shock 
generating  cowl.  The  model  size  and  surface  coordinates  of  the  ramp  and 
cowl  were  obtained  by  scaling  up  by  a  factor  of  2,  Model  E-l  which  was  tested 
earlier  in  the  program  in  the  CAL  and  MSD  shock  tunnel  facilities.  Model  E-3 
ramp,  cowl  contour  geometry  and  overall  dimensions  are  shown  in  the  schematic 
drawing  of  Figure  40. 
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The  Model  E-3  ramp  compression  surface  plate  was  fabricated  from 
6061-T6  Aluminum  in  the  same  manner  as  was  the  Model  C  ramp  surface  plate. 

The  same  approach  to  the  design  of  the  watercooling  system  and  detachable 
tool  steel,  watercooled,  sharp  leading  edge  component  was  used  as  shown  in 
the  photograph  of  Figure  41.  However,  the  Model  E-3  ramp  (and  cowl)  width 
was  reduced  from  is"  (Model  C)  to  14"  to  minimize  tunnel  blockage.  The  NASA 
combination  probe  positioner  plates,  the  .062"  diameter  predrilled  static 
pressure  holes,  and  the  boundary  layer  trip  tapped  holes  are  also  shown  in 
the  figure.  The  Model  E-3  48"  long  ramp  surface  plate  utilized  the  Model  C 
ramp  cast  aluminum  undercarriage  with  its  steel  adapter  section  and  three 
probe  actuators. 

The  Model  E-3  cowl  assembly  consisted  of  a  14"  wide  X  15"  long  6061-T6 
Aluminum  contoured  cowl  plate  (forward  portion  watercooled)  with  a  detachable 
tool  steel,  watercooled,  sharp  leading  edge  component,  and  a  4130  steel  cowl 
support  and  positioning  device  as  documented  in  the  reduced  size  drawing  of 
Figure  42.  The  cowl  positioning  device  was  manually  pre-set  prior  to  each 
test  run  and  provided  for  both  horizontal  and  vertical  cowl  displacement  in 
accomplishing  cowl  angle  of  attack  changes.  Flexible-hose  water  lines 
leading  from  the  cowl  watercooling  system  tubes  provided  the  flexibility 
required  for  cowl  positioning. 

Although  the  Model  C  cowl  was  remotely  positioned  through  angles  of 
attack  in  a  continuous  manner,  a  similar  arrangement  for  the  Model  E  3  cowl 
was  not  practical.  The  Model  C  cowl  was  actuated  through  the  overhead  hatch 
plate  at  test  section  station  43"  which  was  a  considerable  distance  upstream 
(approximately  30")  of  the  shorter  Model  E-3  cowl  leading  edge.  Consequently, 
cowl  drive  shaft  interference  and  stagger  significantly  complicated  the 
efficient  design  of  a  similar  system  for  Model  E-3. 

(2)  Instrumentation 

The  Model  E-3  ramp  was  instrumented  with  80  static  pressure  taps  and 
11  surface-embedded  thermocouples.  The  location  and  identification  of  the 
taps  and  thermocouples  are  documented  in  the  schematic  sketch  of  Figure  43. 
All  static  taps  were  drilled  normal  to  the  curved  surface  of  the  ramp.  The 
method  of  installing  the  .063"  0D  X  ,042  ID  pressure  tubing  and  the  three 
probe  actuators  was  the  same  as  described  for  the  Model  C  ramp.  The 
underside  of  the  Model  E-3  ramp  plate,  however,  was  milled  to  accept  the 
probe  actuators  (at  the  three  probe  stations)  for  normal  to  ramp  curved 
surface  immersion  of  the  NASA-Ames  combination  probe. 

The  Model  E-3  cowl  was  instrumented  with  11  static  pressure  taps  (normal 
to  the  contoured  surface)  and  2  surface-embedded  thermocouples.  The  location 
and  identification  of  the  taps  and  thermocouples  are  documented  in  the 
schematic  of  Figure  44. 

The  three  probe  actuators  and  the  NASA-Ames  boundary  layer  combination 
probe  previously  described  for  the  Model  C  testing  program  were  also  used 
for  Model  e-3  testing.  The  NASA-Ames  combination  probe  reading  stations  for 


Model  E-3  testing  are  indicated  in  Figure  4U  as  are  the  reading  stations 
for  the  “otal  pressure  probe  rakes  A  and  B.  Rake  A  was  a  fixed  rake  while 
Rake  B  was  a  movable  rake,  the  shaft  of  which  attached  to  the  vertical  member 
of  the  cowl  support. 

The  individual  probe  centerline  heights  in  inches  above  the  ramp  surface 
for  each  rake  were  as  follows: 


Rake  A 


Rake  B 


be  No . 

(14  Probes) 

(12  Probes) 

Probe  £  Height 

Probe  (^Height 

1 

.028 

.20 

2 

.070 

.70 

3 

.106 

1.20 

4 

.145 

1.70 

5 

.189 

2.20 

6 

.262 

2.70 

7 

.344 

3.20 

8 

.417 

3.70 

9 

.512 

4.20 

10 

.616 

4.70 

11 

.717 

5.20 

12 

.821 

5.70 

13 

.920 

-- 

14 

1.008 

— 

Data  obtained  with  rakes  A  and  B  were  generally  in  the  field  of  the 
cowl  trailing  edge  expansion  and  thus  were  of  limited  value  in  regard  to 
evaluating  downstream  profile  distortion  effects. 

(3)  Boundary  Layer  Trips 

During  the  Model  E-3  Mach  10  testing  program  a  total  of  16  test  runs 
were  made  with  boundary  layer  trips  installed.  The  boundary  layer  trip 
configuration  selected  for  the  Model  E-3  ramp  is  documented  in  the  schematic 
sketch  of  Figure  45.  Three  rows  of  #2-64  tapped  holes  on  3/8"  centers  were 
installed  across  the  width  of  the  ramp  at  stations  4",  5  1/2",  and  7"  from 
the  ramp  leading  edge  respectively.  Gaps  in  the  spacing  of  the  tapped 
holes  in  rows  2  and  3  were  necessary  to  forestall  breaking  into  the  water 
cooling  chambers  feeding  the  cross-over  water  lines  in  the  aluminum  ramp. 

For  test  runs  3,  4,  10  and  11,  #2-64  stainless  steel  set  screw  trips 
were  installed  in  rows  1  and  2  with  trip  heights  set  at  .040"  and  .060" 
respectively.  For  the  remaining  12  test  runs  utilizing  trips,  the  set 
screw  trips  were  installed  in  rows  1  and  3  with  trip  heights  set  at  .060" 
and  .085"  respectively.  A  typical  installation  of  the  boundary  layer 
trips  (row  1  with  trip  heights  .060"  and  row  3  with  trip  heights  .085") 
is  shown  in  the  photograph  of  Figure  48. 
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(4)  Installation 


A  total  of  35  tost  runs  at  Mach  10.5  and  Mach  7.3  were  made  with  the 
Model  E-3  ramp  assembly  installed  as  shown  in  the  photograph  of  Figure  46. 

The  fixed  axial  orientation  of  the  Model  E  3  ramp  assembly  in  the  test 
section  was  the  same  as  the  axial  orientation  of  the  Model  C  ramp  assembly. 

The  complete  Model  E-3  installation  is  shown  in  the  photograph  of 
Figure  47.  The  cowl  is  shown  with  aluminum  side  plates  installed.  The 
method  of  attaching  the  cowl  support  and  positioning  device  to  the  ramp  and 
cowl  is  documented  in  the  reduced  size  drawing  of  Figure  42.  Cowl  angle  of 
attack  positions  ranging  from  -3°  to  +  18°  were  manually  set  within  .01° 
accuracy  using  a  gunner's  quadrant  graduated  in  mils. 

Test  runs  15  through  41  at  Mach  10.5  and  test  runs  61  through  92  and 
96  through  107  at  Mach  7.3  were  made  with  the  complete  Model  E-3  installation. 
During  test  run  20  a  severe  crack  developed  in  the  ramp  detachable  leading 
edge  component.  The  crack  was  located  near  the  centerline  of  the  model  and 
extended  from  the  sharp  leading  edge  approximately  l"  aft,  then  laterally 
across  the  component  along  a  3,16"  diameter  watercooling  passage.  The  cause 
of  the  failure  was  attributed  to  inadequate  heat  treatment  or  stress  reliev¬ 
ing  of  the  leading  edge  component  during  fabrication.  In  order  to  continue 
testing  with  minimum  delay,  the  damaged  leading  edge  component  was  removed 
from  the  ramp  and  the  Model  C  watercooled  leading  edge  component  was 
installed.  Although  otherwise  interchangeable,  the  Model  C  leading  edge 
component  was  18”  wide  and  provisions  had  to  be  made  as  shown  in  the  photo¬ 
graph  of  Figure  48  to  adapt  to  the  14"  wide  Model  E-3  ramp.  The  transition 
from  the  Model  C  steel  leading  edge  component  surface  to  the  Model  E-3 
aluminum  ramp  surface  was  smooth  and  continuous,  and  the  modified  installa¬ 
tion  was  used  successfully  for  the  remaining  89  test  runs  of  the  program. 

Cowl  side  support  struts,  which  were  installed  early  in  the  cowl-on 
phase  of  test  runs  to  prevent  movement  of  the  cowl  during  tunnel  operation, 
are  shown  in  the  photograph  of  Figure  49.  Not  shown  are  5/b"  wide  stream¬ 
lined  spacer  blocks  installed  between  the  cowl  and  ramp  anu  the  cowl  support 
struts.  These  spacers  were  necessary  to  prevent  local  shock  impingement  from 
the  support  struts  upon  the  ramp,  a  condition  which  could  cause  non-start  or 
choking  of  the  inlet. 

B.  FACILITY  DESCRIPTION 

Models  C,  D  and  E  l  were  tested  in  the  NASA-Ames  3.5  foot  Hypersonic 
Wind  Tunnel  facility  which  is  shown  in  the  sketch  of  Figure  50.  This  tunnel 
is  an  intermittent  blowndown  tunnel  with  running  times  ranging  from  1  to  4 
minutes  and  contoured  nozzle  Mach  number  capabilities  of  5,  7,  and  10.  The 
tunnel  was  operated  at  Reynolds  number  per  foot  varying  from  2  X  10^  to 
3.5  X  10B  at  Mach  number  7  and  5  X  10^  to  2.2  X  106  at  Mach  number  10.  The 
stagnation  pressure  varied  from  50  to  600  psia  at  Mach  number  7  and  from 
400  to  1800  psia  at  Mach  number  10. 
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The  tunnel  operates  with  a  combination  of  initially  high  pressure  air 
and  pre-test  vacuum  pumpdown.  High  pressure  air  from  storage  cylinders  is 
routed  through  the  2000°R  pebble-bed  air  heater  into  the  nozzle  section.  The 
heated  air  is  expanded  through  the  nozzle  (cooled  by  a  helium  injection  system) 
and  enters  the  test  section  of  the  tunnel.  The  air  passes  over  the  model  into 
the  diffuser  section  and  eventually  dumps  into  four  interconnected  large 
vacuum  spheres  which  are  initially  pumped  down  to  l/4  psia  for  tunnel  start¬ 
ing.  The  combination  of  the  high  pressure  air  and  the  vacuum  pumpdown  results 
in  a  pressure  ratio  necessary  for  tunnel  starting.  Model  starting  loads  are 
ht, ] d  to  a  minimum  during  the  tunnel  starting  phase  by  gradually  raising  the 
tunnel  total  pressure.  As  a  safety  precaution  the  tunnel  starting  procedure 
includes  the  sequential  operation  of  a  series  of  interlocks. 

A  view  of  the  tunnel  test  section  with  the  Mach  7  nozzle  section 
installed  is  shown  in  the  photograph  of  Figure  51.  In  the  foreground  of  the 
photograph  is  a  partial  view  of  the  schlieren  system.  The  tl  ree  windows  lo¬ 
cated  between  the  test  section  door  actuator  arms  are  those  u  ed  for  Models 
C,  D  and  E-3  schlieren  photography.  The  downstream  window  frame  contains  the 
special  access  panel  for  watercocling  lines  and  instrumentation  lines.  The 
raised  door  in  the  background  is  the  one  normally  used  for  access  to  the 
test  section  and  model  installation. 

The  tunnel  test  section  is  156"  long  from  the  nozzle  exit  flange  face 
to  the  diffuser  entrance  flange  face.  The  test  section  inlet  diameter  is 
42. ou"  (Station  0")  and  the  exit  diameter  is  45.307"  (Station  156"),  resulting 
in  a  straight  test  section  wall  taper  angle  of  0°  36.5'  for  the(full  156" 
length.  Door  openings  on  each  side  of  the  test  section  are  110  long  (Stat:on 
8"  to  Station  118")  and  taper  in  height  from  21.42"  at  Station  8"  to  22.53" 
at  Station  118".  The  inside  walls  of  the  remotely  operated  doors  are  curved 
and  tapered  to  conform  to  the  interior  walls  of  the  test  section.  The  doors 
also  incorporate  the  schlieren  windows  which  are  centered  at  Stations  24", 

46",  63",  80",  and  102".  The  windows  are  normally  6"  high  and  12"  long,  but 
the  actual  exposed  glass  of  the  windows  measured  5  1 /4 ”  X  11  l/4"  for  this 
testing  program.  Windows  centered  at  Stations  46",  63"  and  80"  (designated 
as  windows  #2,  #3  and  #4  respectively)  were  utilized  for  schlieren  photography. 
A  hatch  opening  6"  wide  X  36"  long  (Station  18"  to  Station  54")  on  top  of  the 
test  section  was  utilized  for  the  special  hatch  plate  for  the  Model  C  cowl 
drive  assembly.  A  facility  hatch  cover,  with  its  interior  surface  conforming 
to  the  interior  wall  curvature  of  the  test  section,  is  installed  for  normal 
testing. 

The  gimbaled  vertical  section  of  the  model  support  assembly  is  vertically 
adjustable  and  enters  the  bottom  of  the  test  section  with  its  centerline  at 
Station  108".  The  8"  0D  horizontal  section  of  the  model  support  assembly 
extends  upstream  to  Station  97.7"  and  accepts  a  keyed  5"  OD  model  strut.  The 
NASA-Ames  Pressure  Transducer  Pod  Housing,  upon  which  Models  C  and  E-3  were 
mounted,  incorporated  a  4.997"  OD  X  2.4"  ID  integral  shaft  with  a  7"  long 
keyed  bearing  surface  which  entered  the  horizontal  model  support  bore. 
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A  single  pass  Z  type  schlieren  system  with  a  mercury  vapor  light  source 
was  used  in  the  Model  C  and  Model  E-3  testing  programs.  The  schlieren  system 
had  provisions  for  simultaneously  viewing  the  model  and  taking  schlieren 
photographs.  Polaroid  photographs  also  could  be  taken  of  the  viewing 
plate  for  immediate  perusal  or  comparison  with  the  test  data. 

Data  acquisition  was  accomplished  by  sampling  of  the  electrical  signals 
from  the  Model  PA208TC  Statham  absolute  pressure  transducers  (lowest  range 
0-5  psia)  on  the  data  tape  of  the  100  channel  high  speed  Beckman  analog  to 
digital  recorder.  Model  surface  pressure  and  temperature  measurements  were 
handled  in  a  standard  manner  by  computing  (using  the  Ajnes  computer  program) 
on  ‘ he  IBM  7o94  the  surface  to  free  stream  static  pressure  ratio,  etc.,  from 
the  beckman  system  data  tape. 

C.  RUN  LOG 

A  run  log  for  tests  conducted  with  Models  C,  D  and  E-3  is  included  hero 
for  reference  purposes.  Test  data  obtained  during  the  Mach  7.3  and  10.5 
testing  at  NASA-Ames  was  compiled  for  use  in  data  analysis  in  References  11 
and  12. 
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RUN  LOC  -  MODEL  E-3  (Cont'd) 
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D.  TEST  RESULTS  AND  ANALYSES 


a.  Model  C 

( 1 )  Flat  Plate  Results 

Initial  data  for  the  model  of  Figure  35  was  obtained  without  the  shock 
generating  cowl  installed,  for  the  purpose  of  determining  the  nature  of  the 
boundary  layer  as  required  for  interpretation  of  subsequent  shock  impingement 
test  results.  The  3  station  boundary  layer  profile  measurements  provided  a 
basis  for  evaluating  regions  of  laminar  and  turbulent  flow,  detailed  boundary 
layer  integral  properties,  and  local  skin  friction  as  discussed  in  the  follow 
ing  sub-sections.  Reduced  boundary  layer  profile  data  are  tabulated  in 
Table  I. 

Wall  Pressure  Data  and  Leading  Edge  Shock  Shape  At  the  lower  values 
of  wall  static  pressure  measured  on  t.ie  3°  wedge  surface  of  Model  C,  it  was 
necessary  to  make  a  small  correction  for  an  apparent  reference  pressure  error 
as  described  in  Section  IV-C.  The  net  result  of  such  a  correction  was  to 
bring  both  the  inviscid  field  ramp  Mach  number  in  line  with  measured  pitot 
data  and  to  obtain  better  agreement  between  measured  wall  static  pressure 
levels  and  values  computed  from  substantiated  theoretical  techniques.  Such 
corrections  had  no  significant  effect  on  higher  values  of  measured  pressure 
such  as  in  the  region  of  a  shock  impingement. 

Using  the  corrected  pressure  levels,  the  leading  edge  shock  shape  was 
computed  using  the  GE/Bertram  technique  of  Reference  1  which  accounts  for 
the  combined  effects  of  bluntness  and  viscous  induced  interactions.  This 
solution  is  superimposed  on  the  schlieren  of  Figure  52  and  the  agreement 
is  seen  to  he  excellent.  Similar  agreement  was  noted  earlier  in  the  program 
for  shock  tunnel  tests  with  Model  A  as  also  reported  in  Reference  1. 

Transition  -  As  a  nominal  assessment  of  the  state  of  the  boundary  layer 
and  of  boundary  layer  trip  effectiveness,  direct  use  was  made  of  the  compress 
ible  profile  power  law  parameter  n  in  the  expression: 

1 

The  flow  was  considered  laminar  for  values  ef  n  from  1  to  2,  transitional  for 
n  from  2  to  5,  and  turbulent  for  values  of  n  greater  than  5.  The  procedure 
for  determining  n  was  to  establish  the  velocity  profile  from  the  probe  data 
using  the  procedure  described  in  Section  IV-C.  The  value  of  n  was  then 
measured  directly  from  the  resulting  velocity  distribution  when  plotted  on 
log-log  paper  as  typified  in  Figures  53  and  54  indicating  that  such  an 
approach  permits  reasonable  analytical  representation  of  the  measured  profile 
data . 
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Summarizing  the  results  indicated  in  Table  I,  at  Mach  10.5,  the  low 
Reynolds  number  data  was  laminar  at  all  three  probe  stations  indicating  that 
natural  transition  would  be  expected  to  occur  at  a  value  of  Re^  greater  than 
950.  The  low  Reynolds  number  data  taken  at  Mach  7.4,  also  indicates  the 
onset  of  boundary  layer  transition  at  values  of  Re^  greater  than  950.  At 
maximum  unit  Reynolds  number,  natural  transition  was  initiated  at  the  first 
and  second  probe  stations  respectively  for  Mach  7.3  and  10.5,  corresponding 
to  values  of  Re Q  of  approximately  1690  and  1360  respectively. 

The  addition  of  boundary  layer  trips  for  high  Reynolds  number  testing 
resulted  in  the  onset  of  transition  ahead  of  the  first  probe  station  at 
Mach  7.3,  but  a  subsequent  decrease  in  the  profile  parameter  n  at  the  second 
and  third  probe  station  suggest  that  the  resulting  boundary  layer  profile 
was  not  completely  in  equilibrium.  On  the  other  hand,  the  boundary  layer 
trips  used  for  high  Reynolds  number  testing  at  Mach  10.5  not  only  promoted 
earlier  boundary  layer  transition,  but  also  permitted  the  attainment  of 
increased  equilibrium  values  of  profile  parameter  n. 

As  a  matter  of  interest,  the  natural  transition  characteristics 
exhibited  in  this  test  data  obtained  with  the  Model  C  are  compared  in  Figure 
55  with  those  determined  from  shock  tunnel  tests  and  heat  transfer  measure¬ 
ments  made  earlier  in  this  program  with  flat  plate  Model  A  and  curved  surface 
Model  E-l  at  Cornell.  (Reference  1.)  The  value  of  n  for  Run  14,  Model  C 
data,  where  the  boundary  layer  profile  is  measured  at  the  third  probe  station 
at  maximum  Reynolds  number,  indicated  the  flow  to  be  turbulent.  At  the  same 
test  conditions,  the  boundary  layer  at  the  second  probe  station  appeared  to 
be  transitional.  Consequently,  in  Figure  55  the  data  point  for  Run  14  is 
plotted  in  such  a  manner  to  indicate  that  boundary  layer  transition  had 
been  completed  at  a  value  of  Re^  of  approximately  2000.  At  this  same  probe 
position,  at  the  lower  unit  Reynolds  number,  the  boundary  layer  was  still 
laminar,  and  thus,  it  is  seen  that  for  both  tests,  transition  occurred  for 
an  Re^  greater  than  1000  and  less  than  2000  on  the  flat  plate  models  as 
compared  to  an  Re^  of  slightly  less  than  1000  on  the  curved  surface  pressure 
gradient  model. 

Variation  of  n  with  Reft  -  To  determine  the  extent  to  which  the  values 
of  n  for  the  hypersonic  non-adiabatic  wall  turbulent  boundary  layer  of  these 
tests  correlate  with  near-adiabatic  wall  results  at  lower  Mach  numbers, 

Figure  56  was  prepared  comparing  the  results  from  Table  1  with  data  presented 
by  Persh  in  Reference  9  and  by  Pinckney  in  Reference  10.  Data  from  this 
program  taken  both  with  and  without  boundary  layer  trips  are  included.  The 
natural  transition  data  not  only  indicate  higher  values  of  n  at  lower  Re,- 
than  the  correlation  of  Reference  9,  but  also  steeper  rate  of  change  of 
n  with  increase  in  Re^.  At  Mach  10.5,  the  addition  of  boundary  layer  trips 
resulted  in  a  significant  increase  in  n  with  only  a  small  change  in  Rep 
from  the  nominal  value  of  2000.  On  the  other  hand,  at  Mach  7.4  with  boundary 
layer  trips  installed,  the  contrasting  decrease  in  n  with  increasing  Re^ 
again  suggests  that  the  boundary  layer  profiles  for  these  particular  test 
conditions  were  not  in  complete  equilibrium.  Consequently,  it  is  implied 
that  the  high  rate  of  mixing  and  turbulence  induced  by  the  boundary  layer 
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trips  at  Mach  7.3  could  not  be  sustained  by  natural  turbulence  levels  exist¬ 
ing  at  these  conditions.  Somewhat  better  correlation  is  obtained  when  n  is 
plotted  as  a  function  of  effective  length  Reynolds  number,  as  also  indicated 
in  Figure  56,  but  the  trend  is  similar. 

Momentum  Thickness  Comparison  -  The  laminar  values  of  momentum  thickness 
as  presented  in  Table  1  have  been  compared  against  the  values  predicted  for 
flat  plate  flow  at  the  measured  boundary  layer  edge  conditions  using  the 
methods  of  Van  Driest  (Reference  13)  and  Cohen-Reshotko  (Reference  M).  T1  is 
comparison  is  shown  in  Figure  57  in  a  manner  also  conducive  to  the  establish 
ment  of  a  laminar  skin  friction  law  based  on  the  test  data.  For  zero 
pressure  gradient  flow,  assuming  a  constant  value  of  n  (as  substantiated  by 
the  profile  measurements)  a  skin  friction  law  can  readily  be  established 
from  the  classical  hypothesis  that  s/x  Cr  (Rex)_z  as  follows: 

-  K  (Re  )~Z  (19) 

x  x 


d0  _  ^f  _  (1-Z)  K 
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From  the  slope  of  a  line  through  the  data  in  Figure  57,  7.  is  determined 
to  be  .50.  K  can  be  evaluated  from  any  point  on  the  curve  to  be  .010.  Thus 
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is  the  experimentally  determined  laminar  skin  friction  law  which  is  seen  to 
be  reasonably  well  predicted  by  the  laminar  theories  of  References  13  and  14. 


A  similar  development  of  a  skin  friction  law  for  the  turbulent  boundary 
laver  data  is  not  Justified  in  that  the  effective  origin  necessarily  varied 
with  test  conditions  and  boundary  layer  trip  configuration.  As  an 
alternative,  use  is  made  of  the  incompressible  skin  friction  law: 
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developed  by  Spence  in  Reference  8  and  adapted  to  compressible  flow  through 
the  Eckert  Reference  Temperature  relation: 
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with  the  result  : 
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where : 
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'Vith  the  aid  of  such  an  expression  and  the  measured  profile  data,  the 
effective  origin  associated  with  each  profile  measurement  was  computed  from: 
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with  the  results  as  reflected  in  Figures  56  and  57.  Equations  21  and  24 
provide  a  Vsis  for  estimating  values  of  laminar  and  turbulent  boundary 
layer  prep*  »  •  l  at  positions  immediately  upstream  of  shock  impingement 
interactic 

Boundary  Layer  Profile  Property  Comparisons  -  For  nearly  all  of  the 
profile  data,  the  compressible  power  law  profile  parameter,  n,  was  calculated. 
In  addition,  for  selected  cases,  the  incompressible  value  of  the  power  law 
profile  parameter,  N,  was  also  calculated.  The  incompressible  cases  were 
calculated  primarily  for  comparison  against  theoretical  profiles.  In 
Figure  58  several  of  the  laminar  profiles  measured  on  Model  C  and  E-3  are 
compared  against  the  Blasius  and  Cohen-Reshotko  flat  plate  solutions.  As 
can  be  seen,  the  experimental  data  falls  below  the  Blasius  solution  indicat¬ 
ing  a  fuller  profile.  Tht re  is  also  a  point  of  inflection  in  the  Model  E-3 
profile  indicating  the  effect  of  the  slight  adverse  pressure  gradient.  It 
should  be  noted,  however,  that  the  experimental  data  can  be  brought  into 
closer  agreement  by  using  a  different  criteria  for  establishing  the  boundary 
layer  edge  conditions.  The  edge  conditions  are  shown  based  on  a  velocity 
ratio  of  .995  of  the  inviscid  wedge  flow  velocity.  Consequently,  with  this 
definition,  any  entropy  layer  due  to  the  leading  edge  shock-boundary  layer 
interaction  is  included  in  the  boundary  layer.  As  an  alternative  approach, 
matching  could  have  been  forced  near  the  wall,  and  boundary  layer  edge 
conditions  (exclusive  of  any  entropy  layer)  defined  as  a  by-product  of  the 
lower  region  profile  matching. 
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Several  comparisons  were  also  made  of  corresponding  compressible  and 


incompressible  profiles  as  indicated  i 
expected,  for  a  typical  turbulent  prof 
profile  parameter  (n)  is  smaller  than 
value  of  n  was  determined  as  indicated 
power  law  profile  parameter  N  for  the 

mined  through  numerical  integration  of 
relation  as  follows: 


and  the  log-log  plotting  procedure  pre 
of  n. 


Figures  59  and  60.  As  would  be 
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iously  described  for  the  determination 


In  Figures  61  to  74,  for  j  irposes  cf  documentation,  typical  boundary 
profile  data  are  plotted  in  terms  of  selected  profile  properties  as  calculated 
in  the  course  of  the  data  reduction. 


From  the  skin  friction  laws  reflected  in  Figure  57,  for  those  cases 
where  the  profile  parameter  n  is  essentially  invariant  with  Reynolds  number, 
the  growth  rate  was  established  for  subsequent  use  in  connection  with 
incident  shock  data  as  follows: 


A  = 


w 

(*  n 

KH 

— 

A 

•J 

(29) 


where  K  and  Z  have  the  previously  es tabl ished  values.  Thus  from  Figure  5/ 
it  is  seen  that  the  hypersonic  laminar  boundary  layer  growth  along  tv  model 
C  ramp  was  in  proportion  to  the  square  root  of  the  Reynolds  number  as  would 
classically  be  predicted  from  low  speed  theory.  For  any  particular  set  of 
similar  turbulent  data  the  growth  rate  can,  of  course,  be  expressed  in  a 
comparable  manner,  but  is  a  function  of  the  particular  value  of  N  associated 
with  that  data. 


Summary  -  From  analysis  of  the  Model  C  ramp  data  it  can  be  concluded 
that  for  Incident  shock  tests  conducted  at  low  te  ;t  facility  unit  Reynolds 
numbers  the  boundary  layer  separation  characteristics  should  be  those  of  a 
laminar  layer  with  consistent  growth  rate  and  local  skin  friction  variations 
derived  from  the  profile  measurements  in  general  agr  nt  with  predicted 
trends.  For  tests  conducted  at  Mach  7.3  and  the  maxima...  facility  unit 
Reynolds  number,  the  incident  shock-boundary  layer  interactions  should  be 
those  associated  with  a  naturally  turbulent  hypersonic  boundary  layer.  At 
other  test  conditions,  boundary  layer  trips  are  required  to  establish 
turbulent  boundary  layer  profiles.  In  that  the  boundary  layer  trips  had 
a  pronounced  effect  on  the  profile  shape,  this  should  be  taken  into 
consideration  in  the  interpretation  of  incident  shock  wave-boundary  layer 
interaction  data.  For  the  hypersonic  boundary  layer  measurements  made 
both  with  and  without  trips,  it  was  found  that  reasonable  representation 
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of  profile  shape  is  possible  through  the  use  of  simple  power  law  profiles. 

(2)  Incident  Shock  Results 

With  the  cowl  installed,  (Figure  35)  interaction  between  the  Model  C 
ramp  boundary  layer  and  the  impinging  cowl  shock  was  studied.  Test  conditions 
and  boundary  layer  trip  geometry  were  predicated  on  earlier  tests  conducted 
without  the  cowl  installed.  At  Mach  10.5,  the  unit  Reynolds  number  was 
regulated  between  approximately  .7  and  2  X  10®  per  foot,  while  tht  cowl 
shock  strength  was  varied  through  the  range  from  0  to  13°  of  turning.  For 
Mach  7.3,  the  unit  Reynolds  number  was  regulated  from  .5  to  3.3  X  10®  and  the 
shock  strength  varied  from  0  to  10°  of  turning.  Nominal  thermocouple  deter¬ 
mined  wall  to  total  temperature  ratios  were  .40  at  Mach  7.3  and  .30  at 
Mach  10.5. 

Flow  Model  Comparisons  -  The  series  of  composite  schlieren-pressure 
distributions  of  Figures  75  to  80  indicate  the  distinct  differences  between 
the  hypersonic  laminar  and  turbulent  boundary  layer  interactions  for 
comparable  incident  shock  strengths.  (Mach  7.3  schlierens  were  selected 
for  these  comparisons  because  of  their  clarity  compared  to  those  obtained 
with  the  lower  density  Mach  10.5  flow.)  The  interaction  region  for  the 
separated  laminar  boundary  layer  is  seen  to  closely  resemble  the  theoretical 
model  of  Section  II-A.  The  shock  reflection  consists  of  a  constant  pressure 
boundary  interaction  followed  by  an  isentropic  compression  back  to  a  flow 
direction  parallel  to  the  wall.  The  isentropic  compression  coalesces  to 
form  the  recompression  shock. 

As  evident  from  the  figures,  the  turbulent  boundary  layer  did  not 
separate  for  the  range  of  shock  strengths  used  in  this  program.  The  turbu¬ 
lent  interaction  region,  in  contrast  to  that  of  the  laminar,  is  characterized 
by  deep  penetration  of  the  boundary  layer  by  the  cowl  shock.  The  downstream 
boundary  layer  is  thinned  by  the  compression  to  a  height  compatible  with 
continuity  and  momentum  constraints.  The  reflected  shock  eminates  from  a 
region  near  the  wall  and  penetrates  through  the  upper  regions  of  the 
boundary  layer  into  the  inviscid  flow  field.  Compared  to  the  laminar  case 
there  is  little  forward  propagation  of  the  resulting  pressure  disturbance 
evident  in  the  wall  static  pressure  data. 

Fr<- m  schlieren  photographs  such  as  those  of  Figures  78  to  80,  measure¬ 
ments  of  the  observable  inviscid  interaction  length  and  boundary  layer  edge 
flow  deflection  were  made  for  comparison  with  the  formulated  flow  model 
discussed  in  Section  II-B.  The  results  as  presented  in  Figures  13  and  81 
indicate  that  the  boundary  edge  flow  deflection  through  the  impinging  shock 
wave,  /'B,  can  be  approximated  by  the  inviscid  shock  deflection  angle  9 
and  the  inviscid  interaction  length  then  estimated  from: 
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based  on  the  flow  geometry. 


Attainment  of  Two-Dimensional  Static  Pressure  Ratio  -  The  overall 
shock  induced  pressure  rise  ratios  are  summarized  in  Figures  82  and  83  for 
the  Mach  7.3  and  10.5  data  respectively.  At  Mach  10.5,  the  Model  C  data  is 
also  compared  against  the  Model  A  shock  tunnel  data  taken  earlier  in  the 
program  and  reported  in  Reference  1.  For  comparison,  corresponding  theo¬ 
retical  pressure  rise  ratios  were  computed  both  inviscidly  and  viscously. 

The  viscous  calculations  were  made  using  the  G.E. /Bertram  solution  to  obtain 
the  effective  shock  strength  due  to  the  viscous  and  bluntness  induced  lead¬ 
ing  edge  Interaction.  The  effective  shock  strength  was  then  used  to  calculate 
the  final  pressure  rise.  The  data  agrees  very  well  with  the  estimates; 
generally  falling  between  the  limits  set  by  the  invlscid  and  viscous  solutions, 
both  for  the  separated  and  attached  flow  cases. 

Laminar  Theory/Data  Separation  Plateau  Comparisons  -  The  hypersonic 
laminar  separation  data  was  compared  tc  the  Erdos  and  Pallone  formulations, 
which  give  good  agreement  at  lower  Mach  numbers,  for  the  prediction  of 
plateau  levels  and  interaction  lengths  (Reference  7).  For  this  purpose, 
the  series  of  pressure  distributions  shown  In  Figures  84  to  93  were  prepared. 

At  both  Mach  7.3  and  10.5,  separation  plateaus  are  seen  to  form  at  cowl 
angles  of  3  to  4°  (Figures  84,  85  and  90).  As  the  shock  strength  is 
increased,  the  length  of  the  pressure  plateau  increases.  Superimposed  on 
these  pressure  distributions  is  the  pressure  plateau  levels  predicted  by 
the  theory  of  Erdos  and  Pallone  using  the  laminar  skin  triction  law 
established  from  cowl-off  testing  of  this  model  as  reported  in  the  previous 
section.  The  excellent  agreement  is  taken  to  indicate  that  the  Erdos  and 
Pallone  semi-empirical  formulation: 


based  on  low  Mach  number  data,  accurately  predicts  laminar  separation 
plateau  pressure  levels  at  hypersonic  speeds. 

Laminar  Separation  Length  Theory/Data  Comparisons  -  In  Reference  7 
Erdos  and  Pallone  formulated  an  expression  for  the  separation  length  which 
correlates  well  with  low  Mach  number  data  in  terms  of  boundary  layer  thick¬ 
ness  and  final  pressure  rise  minus  plateau  pressure  rise.  In  comparing 
this  formulation  with  the  present  data,  a  strong  Mach  number  dependence 
was  found  as  indicated  in  Figure  94.  To  further  evaluate  the  validity  of 
this  apparent  Mach  number  dependence,  two  recent  sets  of  data  from 
independent  sources  were  examined.  As  also  indicated  in  Figure  94,  the 
Reference  15  data,  at  a  nominal  Mach  number  of  10  and  at  a  length  Reynolds 
number  of  .5  X  10®  (compared  to  1.5  X  10®  for  the  present  Mach  10.5  test) 
agree  with  the  trend  of  the  data  obtained  in  this  program.  Reference  3 
data,  on  the  other  hand,  generally  is  in  agreement  with  the  Erdos  and 
Pallone  formulation.  It  was,  however,  taken  at  Mach  2.21,  closely  corres¬ 
ponding  to  the  conditions  of  the  original  Erdos  and  Pallone  correlation. 
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Failure  to  observe  the  apparent  Mach  number  dependence  evident  in  Figure  94 
when  predicting  laminar  separation  lengths  is  seen  to  result  in  substantial 
over-estimation  at  hypersonic  speeds.  Admittedly,  any  shock  induced  boundary 
layer  transition  effects  present  in  the  data  could  be  a  factor  in  the  reduced 
separation  lengths,  but  this  appears  somewhat  unlikely  considering  the 
consistent  pattern  of  the  present  data  and  the  agreement  with  the  experiment 
of  Reference  15. 

Laminar  Theory/Data  Pressure  Distribution  Comparison  -  Results  obtained 
with  the  present  Laminar-Shock/Boundary  Layer  Interaction  program  discussed 
in  Section  II-A  were  compared  with  the  Model  C  laminar  interaction  data  at 
both  Mach  10.5  and  7.3.  Since  profile  measurements  indicated  none  of  the 
separated  data  remained  laminar  downstream  of  the  reattachment  point,  only 
two  runs  with  relatively  weak  interaction  give  evidence  of  wholly  laminar 
flow  and  could  be  compared  directly  with  the  theory.  Run  61  at  Mach  7.3  was 
the  only  run  with  a  laminar  boundary  layer  measured  downstream  of  the  inter¬ 
action  (probe  station  3  was  approximately  17  inches  downstream  of  the  shock 
impingement  point)  and  Run  50  at  Mach  10.5  had  a  laminar-transitional 
interaction.  Data  from  both  of  these  runs  are  compared  with  the  computer 
program  calculations  in  Figures  95  and  96.  The  predicted  forward  propagation 
of  the  pressure  and  the  initial  slope  of  the  pressure  rise  appear  to  be  in 
very  good  agreement  with  the  experimental  data.  The  overall  pressure  rise  is 
matched  near  the  end  of  the  plate,  however,  there  is  a  definite  overshoot  of 
the  experimental  pressure  data  in  the  interaction  region  which  is  not 
theoretically  predicted.  (The  maximum  pressure  level  in  this  region  could 
have  been  more  nearly  approximated  theoretically  by  using  an  increased  shock 
strength  such  as  might  be  associated  with  a  wave  generated  by  a  viscous  cowl. 
A  correction  of  less  than  .4  degrees  to  the  shock  induced  turning  angle 
would  increase  the  final  theoretical  pressure  rise  until  it  matched  the 
experimental  data  peak.) 

Laminar  Theory/Data  Mass  Flow  Comparisons  -  Figures  97  and  98  reflect 
the  theoretical  mass  flow  variation  through  the  interaction  region  for  the 
laminar-laminar  case  (Run  61)  and  the  laminar-transitional  case  (Run  50) 
as  determined  by  the  Laminar-Shock/Boundary  Layer  Interaction  program 
(laminar-laminar  calculation).  The  parameter  plotted  in  Figures  97  and  98 
is  a  mass  flow  parameter  related  to  the  actual  mass  flow  by  the  following 
expression : 


m  "  Wi 


(32) 


Where:  m  =  mass  flow  in  boundary  layer  slug/sec 

W__  =  mass  flow  parameter  -  from  Laminar-Shock/ 
Boundary  Layer  Interaction  program 


BL 


=  density  at  boundary  layer  edge  before 
Interaction 

u =  velocity  at  point  (1) 

p  =  local  density  at  boundary  layer  edge 
e 

ug  =  local  velocity  at  boundary  layer  edge 
6  =  local  boundary  layer  thickness 

6*  =  local  displacement  thickness 


Figures  97  and  98  indicate  a  slight,  temporary  decrease  in  boundary 
layer  mass  during  the  interaction.  This  could  be  caused  by  the  similar 
solution  approximation  used  in  the  formulation  of  the  Laminar-Shock/Boundary 
Layer  Interaction  program.  The  decrease  in  mass  occurs  in  the  free  inter¬ 
action  region  and  is  most  apparent  for  the  incipient  separation  cases.  This 
discrepancy  does  not  appear  to  have  a  significant  effect  on  the  overall 
mass  flow  predictions,  however. 

The  theoretical  ratio  of  the  mass  flow  downstream  to  the  mass  flow 
upstream,  CA^/ A2 )  was  compared  with  the  boundary  layer  mass  flow  ratio 
measured  across  the  interaction  region  during  the  experimental  program. 

For  the  laminar-laminar  interaction  the  mass  flow  ratio  across  the  inter¬ 
action  agrees  very  well  with  the  ratio  predicted  theore.  ally.  From  the 
experimental  data  obtained  during  Run  61  the  mass  flow  ratio  across  the 
Interaction,  measured  at  probe  stations  number  1  and  2  was  determined  to 
be  1.28  whereas  a  theoretical  value  of  1.294  at  the  same  stations  is 
predicted  using  the  Laminar-Shock/Boundary  Layer  Interaction  program. 

For  the  laminar-transitional  interaction  (Run  50) ,  the  experimentally 
determined  mass  flow  ratio  across  the  interaction  was  much  higher,  confirm¬ 
ing  that  this  interaction  was  transitional  (2.95  experimentally  as  compared 
to  1.71  theoretically). 

Laminar  Theory/Data  Displacement  Thickness  Comparison  -  A  theory/data 
comparison  of  the  displacement  thickness  at  the  aft  probing  station  indicated 
very  favorable  agreement  between  the  shock  boundary  layer  prediction  and  the 
reduced  experimental  data.  For  the  laminar-laminar  case,  the  displacement 
thicknesses  for  the  theoretical  and  experimental  cases  are  .134  and  .138 
inches  respectively.  For  the  laminar-transitional  case,  the  displacement 
thicknesses  are  .190  and  .220  inches  respectively. 

Theoretical  Interaction  Region  Laminar  Skin  Friction  -  In  Figures  99 
and  100  the  shear  parameter  (£)  is  plotted  versus  axial  location  for  the 
Mach  10.5  and  7.3  test  conditions  used  in  the  previous  comparisons.  The 
shear  parameter  is  defined  as  follows: 
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velocity  gradient  at  wall 


Re 

9  =  Reynolds  number  based  on  0 

Cf  z  skin  friction  coefficient  based  on  boundary 
layer  edge  conditions 

\  -  viscosity  parameter  z 


Separation  is  indicated  by  negative  values  of  l. 


As  can  be  seen  from  these  plots,  incident  shock/boundary  layer  separa¬ 
tion  was  not  predicted  in  either  case,  as  indicated  by  the  positive 
value  of  the  skin  friction  parameter  throughout  the  interaction  region.  For 
all  of  the  separated  flow  cases,  test  data  analysis  indicates  that  the 
boundary  layer  had  transitioned  downstream  and  was  fully  developed  turbulent 
flow  by  the  aft  probing  station.  Consequently,  theory/data  comparisons  for 
those  conditions  are  not  applicable. 


Turbulent  Theory /Data  Separation  Comparison  -  With  upstream  turbulent 
boundary  layer,  no  observable  separation  was  noted  either  in  the  pressure 
distributions  or  in  the  schlieren  pictures.  Based  upon  pre-test  estimates, 
turbulent,  free  interaction  pressure  plateau  levels  of  the  order  of: 


(34) 


were  anticipated,  based  on  estimates  made  using  the  method  of  Reference  7. 
Theoretical  plateau  pressure  ratios  (based  upon  the  turbulent  flow  conditions 
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?<p  =  1800  psla,  Tt  z  1950°R,  M  =  10.5)  of  approximately  4.3  were  calculated. 
This  comparts  to  the  maximum  final  pressure  rise  ratio  of  40  measured  during 
Run  53  when  the  cowl  was  set  at  a  maximum  pngle  of  13°  and  no  separation 
plateau  was  observed.  More  detailed  analysis  subsequently  delineated  the 
importance  of  upstream  shape  factor  on  such  separation  predictions  as  reflected 
in  the  theoretical  curves  of  Figure  19  in  Section  II-B.  The  turbulent  shock/ 
boundary  layer  flow  model  developed  in  Section  II-B  was  used  for  post-test 
estimates  of  the  onset  of  boundary  layer  separation.  This  procedure  indicated 
that  the  boundary  layer  should  not  be  expected  to  separate  at  the  Model  C  test 
conditions.  Limited  additional  theory/data  comparisons  were  made  using  test 
data  from  the  literature,  wherein  both  upstream  profile  data  and  wall  static 
pressure  data  (from  which  the  nature  of  the  flow;  separated  versus  attached 
could  be  assessed)  were  reported.  The  results  are  summarized  in  Figure  101, 
suggesting  the  method  of  Section  II-B  is  an  acceptable  engineering  approxima¬ 
tion  for  assessing  imminent  boundary  layer  separation.  In  that  the  computed 
results  are  sensitive  to  upstream  profile  shape  and  wall  to  total  temperature 
ratio,  it  is  important  that  reasonably  accurate  measurements  of  the  parameters 
be  made  as  indicated  in  the  parametric  results  of  Figure  102. 

Boundary  Layer  Thickness  Ratio  Theory/Data  Comparison  -  Results  obtained 
for  estimates  of  the  boundary  layer  thickness  change  across  the  shock  impinge¬ 
ment,  using  the  turbulent  boundary  layer  formulation  presented  in  Section  II-B, 
are  compared  against  data  taken  during  the  Model  C  Mach  7.3  tests.  Measure¬ 
ments  of  boundary  layer  thickness  on  either  side  of  the  interaction  were  made 
from  schlieren  pictures  such  as  those  in  Figures  78  to  80.  As  can  be  seen 
from  Figure  103,  there  is  good  agreement  between  the  theory  and  the  data. 

The  importance  of  properly  accounting  for  the  shock  induced  downstream  profile 
distortion  is  also  apparent  from  the  comparison  of  Figure  103.  The  assumption 
of  a  constant  profile  parameter  n  across  the  immediate  interaction  region  is 
seen  to  result  in  a  significant  over-estimation  of  downstream  profile  thinning 
at  the  higher  shock  strengths.  The  relative  agreement  obtained  when  profile 
distortion  is  accounted  for  (Reference  5)  is  considered  indirect  substantia¬ 
tion  for  the  supposition  that  mass  flux  rates  measured  downstream  of  the 
immediate  interaction  did  not  occur  across  the  interaction,  but  are  rather 
the  result  of  subsequent  downstream  mixing. 

Turbulent  Interaction  Mass  Flux  -  The  boundary  layer  profile  measurements 
and  integral  properties  established  from  cowl-on  Model  C  testing  are  summarized 
in  Table  II.  Since  the  total  test  conditions  varied  slightly  between  wind 
tunnel  runs,  the  mass  flow  was  subsequently  normalized  to  a  standard  set  of 

test  conditions.  The  normal  total  conditions  for  Mach  10.5  are  a  total 

pressure  of  600  psia  and  a  total  temperature  of  1800°R  for  the  low  Reynolds 
number  runs,  and  1800  psia  and  2000°R  respectively  for  the  high  Reynolds 

number  runs.  At  Mach  7.3,  the  normal  total  conditions  are  100  psia  and 

1400°R  for  the  low  Reynolds  number  and  600  psia  and  1500°R  for  the  high 
Reynolds  number  runs.  The  following  expressions  were  used  to  determine  the 
normalized  mass  flow: 

m  =  m  Laminar  Flow  (35) 

n  Tn'  t 


4U 


and 


An  =  m  ^TT/T-»n^*6  ^PTn/PT^'8  Turbulent  Flow  (36) 

Where:  m  =  normalized  mass  flow 

n 

• 

m  =  measured  mass  flow 

T  =  normalized  total  temperature 

Tt  =  total  temperature 

=  normalized  total  pressure 
Tn 

P^  =  total  pressure 

The  above  two  expressions  are  based  upon  the  assumption  that  the  boundary 
layer  grows  as  the  x*®  and  xf®  for  the  laminar  and  turbulent  case  respectively; 
where  x  is  the  axial  distance  along  the  plate.  The  cowl-off  data  analysis  of 
the  preceding  section  bears  out  the  laminar  growth  rate,  whereas  considera¬ 
tion  of  effective  origin  and  variation  on  n  reduce  the  accuracy  of  such  an 
assumption  for  the  turbulent  data. 

Profile  measurements  made  downstream  of  the  shock  impingement  point, 
at  probe  station  3,  with  upstream  turbulent  flow  also  indicate  significant 
mass  flux  occurred  compared  to  the  normal  boundary  layer  growth  rate  upstream 
of  the  shock  impingement  point.  This  is  summarized  in  Figure  104  indicating 
a  general  increase  in  mass  flux  as  shock  strength  is  increased.  Comparative 
calculations  were  made  using  the  downstream  boundary  layer  edge  conditions 
to  establish  a  new  unit  mass  flow  and  a  new  effective  origin  for  the  down¬ 
stream  flow  based  on  the  boundary  layer  thickness  immediately  downstream 
of  the  shock  interaction  (using  the  technique  of  Section  II-B)  and  assuming 
a  downstream  growth  rate  at  the  new  effective  Reynolds  number  similar  to 
that  for  the  upstream  profile.  Such  calculations  were  found  to  consistently 
underestimate  the  mass  flux  as  indicated  by  the  station  3  probe  measure¬ 
ments  and  are  taken  as  an  indication  that  the  downstream  profile  distortion 
should  be  considered  in  more  detailed  analysis  of  downstream  mass  flux.  It 
should  be  noted,  however,  that  the  Mach  7.3  data  point  (Run  58)  does  not 
outwardly  appear  to  agree  with  the  Mach  10.5  curve.  This  apparent  discrepancy 
may  be  caused  by  the  fact  that  for  the  Mach  10.5  case  the  probing  station  was 
9"  downstream  of  the  interaction  region  compared  to  approximately  2"  for 
Mach  7.3,  suggesting  that  the  greatest  mass  flux  occurred  downstream  of  the 
immediate  interaction  region.  This  is  in  direct  contrast  to  the  observation 
that  the  Mach  7.3  and  Mach  10.5  laminar  data  do  agree  reasonably  well  even 
though  there  is  the  sa ~,e  substantial  difference  between  the  location  of  the 
probing  station  with  respect  to  the  interaction  region,  suggesting  that  the 
transition  induced  mass  flux  occurred  in  the  immediate  region  of,  or  shortly 
downstream  of  the  interaction  region. 
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Profile  Change  Across  the  Shock  Interaction  -  A  comparison  of  the  shock 
interaction  profile  data  with  the  flat  plate  data  as  typified  in  Figure  105 
Indicates  the  nature  of  the  profile  change  across  the  shock  interaction. 

For  the  turbulent  interaction  at  Mach  10.5,  the  downstream  value  of  n,  at 
the  probe  station  3,  is  approximately  the  same  as  that  of  the  upstream 
profile.  The  corresponding  shape  factor  decrease  is  thus  primarily  due  to 
a  decrease  in  Mach  number  across  the  interaction.  This  is  illustrated  from 
the  data  of  Runs  19  and  48.  The  experimental  ratio  of  the  shape  factors 
<e19/h48>  1®  1*89*  Assuming  a  constant  profile  parameter,  n,  using  the 
measured  variation  in  boundary  layer  Mach  number,  the  theoretical  ratio  is 
1.96.  It  should  again  be  noted  that  the  probing  station  at  Mach  It). 5  is 
approximately  9"  downstream  from  the  interaction  region  which  explains  why 
a  disturbed  profile  was  not  measured.  For  the  Mach  7.3  case,  however,  the 
probing  station  was  closer  to  the  interaction  region  (approximately  2" 
downstream).  A  plot  of  the  flat  plate  and  incident  shock  profiles  for  this 
condition  is  shown  in  Figure  106  and  is  evidence  that  the  profile  has 
already  adjusted  into  a  higher  power  law  type  profile  instead  of  a  more 
distorted  profile  requiring  more  than  a  single  parameter,  such  as  n,  for 
reasonable  analytical  representation.  Consequently,  such  data  provides 
additional  substantiation  for  the  turbulent  flow  formulation  of  Section  II-B 
and  of  Reference  5. 

To  graphically  illustrate  and  document  the  results  of  Table  II,  a 
series  of  summary  figures  were  prepared  (Figures  107  to  113)  which  show  the 
velocity  profile  both  upstream  and  downstream  of  the  shock  interaction  with 
respect  to  the  pressure  distribution.  The  mass  flow,  boundary  layer  edge 
Mach  number,  and  velocity  profile  parameter  are  also  indicated  in  these 
figures.  (The  velocity  profiles  presented  have  been  plotted  from  computer 
program  data  reduction  results  based  on  actual  profile  measurements.) 

Summary  -  From  -the  analysis  of  the  incident  shock  Model  C  data  it  can 
be  concluded  that  the  hypersonic  laminar  boundary  layer  separation  pressure 
plateau  levels  can  be  accurately  predicted  by  the  method  of  Erdos  and  Pallone, 
whereas  the  measured  interaction  lengths  exhibit  a  significant  Mach  number 
dependency  which  must  be  accounted  for.  Laminar  separation  was  found  to 
occur  at  a  shock  induced  pressure  ratio  of  5  at  Mach  10.5.  The  apparent 
stability  of  the  hypersonic  turbulent  boundary  layer,  as  documented  by  the 
attached  flow  characteristics  in  spite  of  shock  pressure  rise  ratios  of  40 
at  Mach  10.5  and  13  at  Mach  7.3  is  in  accord  with  the  formulation  of  Section 
II-B  but  was  not  anticipated  on  the  basis  of  lower  Mach  number  correlations. 
Limited  comparisons  with  turbulent  boundary  layer  separation  data  obtained 
by  other  investigators  attest  to  the  engineering  usefulness  of  the  formula¬ 
tion  of  Section  II-B  for  predicting  imminent  separation.  It  can  thus  be 
concluded  that  analyses  of  the  incident  shock/turbulent  flow  data  have 
provided  the  sought  basis  for  a  formulation  of  an  appropriate  flow  model. 

b.  Model  D 

Flow  Model  Formulation  -  The  purpose  of  this  test  was  to  examine  the 
effect  of  swept  vane  generated  shocks  on  the  Model  C  r3mp  laminar  and 
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turbulent  boundary  layers.  Measurements  consisted  of  ramp  centerline  and 
vane  surface  pressures,  and  ramp  pitot  profile  data  taken  at  Model  C  probe 
stations  2  and  3  behind  the  swept  shock  "x"  wave  intersection.  Data  were 
obtained  for  vane  shock  angles  of  6  and  8°.  A  simplified  flow  model  was 
formulated  for  analysis  of  the  inviscid  compression  system. 

For  swept  wedges  of  infinite  span,  that  is  without  tip  effects,  a  simple 
but  accurate  estimation  of  the  chordwise  pressure  distribution  including 
viscous  and  bluntness  effects  can  be  made  based  on  the  G.E. /Bertram  technique 
discussed  in  Reference  1.  From  simple  inviscid  sweep  theory,  it  can  be 
shown  that  the  deviation  in  direction  with  respect  to  the  free  stream,  of  the 
resultant  ramp  velocity  vector  is  expressed  by: 


o 


Calculations  indicate  that  for  most  hypersonic  inlet  applications,  this  angle 
is  small  (of  the  order  of  several  degrees)  and  thus  suggest  that  strip  theory 
should  be  adequate  in  the  absence  of  significant  cross  flow  effects.  Conse¬ 
quently,  for  the  case  of  sharp  leading  edges,  the  G.E. /Bertram  analysis  can 
be  utilized  directly.  For  the  more  practical  case  of  blunted  leading  edges, 
the  pressure  increment  induced  by  bluntness  can  again  be  estimated  using  the 
G.E. /Bertram  method,  but  with  the  leading  edge  drag  coefficient  corrected 
for  sweep  effects: 

C_  =  C_  cos2(i  (38) 

dg  Dfi  =  o 

as  is  also  suggested  by  simple  sweep  theory.  The  same  iterative  calculation 
procedure  developed  for  the  unswept  leading  edge  with  combined  bluntness  and 
viscous  interaction  can  then  be  applied. 

For  the  case  of  a  finite  span  swept  ramp  such  as  that  represented  by 
vertical  installation  on  a  compression  ramp  such  as  Model  D,  an  additional 
boundary  condition  must  be  considered  to  treat  the  flow  influenced  by  the 
finite  span.  For  this  case,  inviscid  analysis  based  on  simple  sweep  theory 
can  be  modified  to  assure  that  the  flow  vector  downstream  of  the  swept 
vertical  ramp  compression  surface,  in  the  vicinity  of  the  corner,  is  not 
only  in  the  plane  of  the  vane,  but  is  also  in  the  plane  of  the  horizontal 
compression  surface  as  well.  The  flow  model  shown  in  Figures  114  and  115 
was  hypothesized  for  establishing  inviscid  flow  field  properties  such  as 
required  in  establishing  boundary  layer  edge  conditions.  This  formulation 
represents  an  extension  of  inviscid  swept  wedge  theory  in  that  it  also 
handles  associated  corner  flow,  as  is  found  in  Model  D.  The  flow  direction 
boundary  conditions  for  the  Model  D  swept  wedge-corner  combination  can  be 
satisfied  by  hypothesizing  the  existence  of  a  downstream  two-dimensional 
centered  expansion  fan  originating  along  the  junction  line  between  the  vane 
generated  shock  sheet  and  the  ramp  surface.  With  such  an  expansion  fan  the 
final  flow  vector  can  be  positioned  both  :in  the  plane  of  the  vane  and  the 
ramp.  Figure  114  depicts  this  "shock-expansion"  flow  model  and  lists  the 
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equations  needed  to  calculate  this  flow  angle. 

Vane  and  Main  Ramp  Pressure  Data/Theory  Correlation  -  The  measured  vane 
pressure  distributions  shown  in  Figure  116  support  the  flow  model  of  Figures 
114  and  115.  Pressure  levels  upstream  of  the  expansion  fan  are  predicted 
by  strip  theory  and  downstream  pressure  levels  tend  toward  the  values 
associated  with  the  shock-expansion  approach.  In  the  intermediate  region, 
preasures  reach  a  common  level  by  communication  via  the  subsonic  flow  in  the 
boundary  layer.  Of  interest  is  the  "lag"  of  the  upper  row  pressures  in  all 
Instances.  Since  both  upper  and  lower  sets  of  orifice  locations  have  been 
normalized  with  respect  to  their  distances  from  the  vane  leading  edge  (see 
symbol  key  insert  on  Figure  116)  the  "lag"  represents  evidence  of  the  physical 
existence  of  a  downstream  expansion. 

Measured  ramp  centerline  pressure  distributions  are  shown  in  Figures  117, 
118  and  119.  Superimposed  on  the  data  are  the  invlscid  theoretical  "x"  wave 
pressure  levels.  Again  go'  agreement  with  the  modified  shock  expansion 
theory  is  obtained.  Strip  «.  .eory  (about  the  same  pressure  level  as  swept 
wedge  theory)  over-estimates  the  measured  pressure  level  to  such  an  extent 
that  it  becomes  inadequate  for  handling  the  swept  wedge-corner  flow  problem. 

Of  interest  is  the  three-dimensionality  of  the  vane  shock  generated  flow 

field  which  is  Indicated  to  some  extent  by  the  relative  levels  of  the  off- 

centerline  ramp  pressures.  In  using  the  shock-expansion  technique  a  correc¬ 
tion  in  the  downstream  displacement  thickness  slope  of  1  for  both  6°  vane 
settings  and  1.5°  for  the  8°  vane  setting  was  made.  These  values  represent 
nominal  estimates  based  on  ramp  boundary  layer  profile  measurements  made  both 
with  and  without  vanes. 

As  can  be  seen  from  the  preceding  figures,  no  vane  shock  induced  boundary 
layer  separation  was  observed  in  the  case  of  the  6°  and  8o  vane-turbulent 
flow  tests.  For  the  6°  vane  setting,  the  laminar  boundary  layer  did  separate 
as  evidenced  in  Figure  119.  As  in  the  case  of  the  Model  C  data,  the  Erdos 
and  Pallone  theory  predicts  the  correct  laminar  plateau  level.  The  separa¬ 
tion  region  in  this  case  is  longer  than  that  observed  on  Model  C  for  a 

similar  downstream  pressure  rise.  The  complex  cross  flow  induced  three- 
dimensionality  of  this  problem  necessitates  the  acquisition  of  data  at  other 
conditions  before  a  comprehensive  correlation  with  two-dimensional  flow 
separation  lengths,  such  as  determined  in  Model  C  testing,  can  be  obtained. 

Main  Ramp  Profile  Data  Analysis  -  Profile  data  shown  in  Figures  120, 

121,  122  and  123  again  demonstrate  the  complexity  of  the  vane  shock  "x"  wave 
flow  field.  Since  the  static  probe  measurements  were  of  limited  value,  only 
first  order  corrections  for  normal  pressure  gradients  (such  as  when  the 
probe  passed  through  a  vane  shock)  could  be  made  in  reducing  the  pitot  probe 
data  to  Mach  number  profiles.  The  validity  of  the  static  pressure  corrections 
is  supported  by  the  agreement  between  measured  boundary  layer  edge  Mach 
numbers  and  predictions  based  on  the  shock-expansion  theory.  Increases  in 
boundary  layer  thickness  and  mass  flow  through  the  vane  shock  system  are 
shown  in  Figure  124. 
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The  profile  data  taken  in  the  vane  shock  field  indicate  that  the 
rate  of  mass  flux  into  the  boundary  layer  increases  as  the  vane  angle 
increases.  This  growth  rate,  in  excess  of  the  normal  flat  plate  rate,  is 
attributed  to  the  reduction  in  downstream  ramp  wetted  area  as  well  as  the 
reduced  effective  Reynolds  number  across  the  vane  shock  system.  Downstream 
of  the  vane  generated  shock  system  local  lower  Mach  number-higher  pressure 
areas  tend  to  counteract  the  boundary  layer  thickening  associated  with  this 
increased  axial  mass  flux.  The  overall  boundary  layer  thickening  observed 
through  this  swept  vane  compression  system  is  contrary  to  the  thinning 
generally  observed  in  tests  conducted  with  the  two-dimensional  incident 
shock  Model  C.  Future  formulation  of  techniques  for  analyzing  the  boundary 
layer  characteristics  should  provide  for  mass  addition  across  the  shock 
system  resulting  from  mixing  rates  associated  with  profile  distortion  as 
well  as  that  due  to  effective  origin  shift  and  change  in  boundary  layer 
edge  unit  Reynolds  number. 

Summary  -  In  evaluating  the  effect  of  the  Model  D  swept  vane  shock 
induced  pressure  rise  on  the  Model  C  compression  ramp  boundary  layer,  it  was 
found  necessary  to  extend  swept  wedge  flow  theory  to  also  handle  the  Model  D 
type  corner  flow.  This  modification  in  the  theory,  here  called  the  shock- 
expansion  approach,  yields  results  well  substantiated  experimentally  by  vane 
and  ramp  surface  pressure  levels. 

For  the  6°  and  8°  vane  settings  tested  no  turbulent  flow  separation  was 
observed.  This  result  is  in  agreement  with  the  separation  theory  formulated 
elsewhere  in  this  report,  which  shows  that  the  generated  pressure  rise, 
lower  than  obtained  for  Model  C,  is  not  great  enough  to  induce  boundary 
layer  separation  for  the  relatively  high  profile  parameter  cold  wall  boundary 
layer  generated  on  the  Model  C  ramp. 

Laminar  separation  was  achieved  with  the  6°  vane  setting.  Erdos  and 
Pallone  theory  predicts  the  correct  plateau  pressure  level.  The  separation 
region  in  this  case  is  longer  due  to  cross  flow  than  that  observed  on  Model 
C  for  a  similar  downstream  pressure  rise. 

Profile  data  taken  downstream  of  the  swept  vane  shock  "x"  wave  location 
indicate  that  future  formulation  of  techniques  for  analyzing  inlet  compression 
systems  of  this  type  should  provide  for  mass  addition  across  the  shock  system 
resulting  from  mixing  rates  associated  with  profile  distortion  as  well  as 
that  due  to  effective  origin  shift  and  change  in  boundary  layer  edge  unit 
Reynolds  number. 

c.  Model  E-3 

(1)  No-Cowl  Configuration 

The  Model  E-3  cowl-off  test  data  was  analyzed  in  the  same  manner  as 
that  of  Model  C  with  the  following  exception.  As  pointed  out  in  Section 
IV-C  variable  static  pressure  normal  to  the  ramp  surface  was  used  in  reduc¬ 
ing  the  boundary  layer  profile  data  compared  to  the  use  of  constant  static 
pressure  in  the  Model  C  data  reduction.  In  Figures  125  to  127  the  measured 
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static  pressure  variations  normal  to  the  wall  at  each  of  the  three  probe 
stations  is  compared  with  a  Method  of  Characteristics  solution  obtained  using 
the  lnviscid  ramp  coordinates.  As  can  be  seen  from  these  figures,  a  signifi¬ 
cant  normal  pressure  gradient  existed  at  the  aft  probing  station.  The 
theoretical  and  measured  gradients  at  all  three  probe  stations  are  in  reason¬ 
able  agreement,  suggesting  the  applicability  of  such  an  analytical  approach 
in  the  absence  of  measured  normal  pressure  gradients. 

Pressure  Distributions  -  A  typical  measured  ramp  pressure  distribution 
fc»  the  cowl-off  configuration  at  Mach  10.5  is  compared  to  the  pre-test 
estimates  and  also  to  the  post-test  Method  of  Characteristics  solution  in 
Figure  128.  As  can  be  seen  from  the  figure,  excellent  agreement  is  obtained 
in  the  three-way  comparison.  A  slight  waviness  in  the  model  wall,  consistent 
with  manufacturing  tolerances,  is  evident  in  the  pressure  data  and  also  in 
the  Method  of  Characteristics  solution  which  is  based  on  the  actual  model 
coordinates. 

Boundary  Layer  Integral  Properties  -  Results  of  Model  E-3  profile  data 
reduction  are  tabulated  in  Table  III  for  both  the  Mach  10.5  and  7.3  test 
conditions,  including  data  obtained  wi?h  i  ,e  cowl  installed.  The  latter 
represent  profile  measurements  taken  oni>  upstream  of  the  cowl  shock  impinge¬ 
ment  point  in  that  model  internal  geometry  constraints  precluded  measurements 
of  flow  profiles  downstream  of  the  shock  interaction.  A  series  of  comparisons 
of  these  boundary  layer  integral  properties  are  shown  in  Figures  129  to  131. 

In  Figure  129,  experimentally  determined  momentum  thickness  versus  Re  is 
compared  against  the  theories  of  Van  Driest  and  Cohen  and  Reshotko  for  flat 
plate  skin  friction.  As  would  be  expected  on  the  basis  of  the  Model  C 
results,  there  is  a  very  favorable  comparison  between  the  theories  and  the 
laminar  boundary  layer  data.  The  turbulent  data,  however,  deviates  from 
the  theoretical  curve  as  would  be  expected  because  of  the  difference  between 
turbulent  and  laminar  skin  friction.  An  interesting  feature  of  this 
comparison  is  that  the  laminar  flat  plate  theories  agree  well  with  the 
experimental  data  obtained  in  the  mild  adverse  pressure  gradient  field  of 
Model  E-3.  It  should  be  noted,  however,  that  Re  was  calculated  from  the 
local  Reynolds  number  and  the  length  to  the  leading  edge  of  the  plate. 
Actually,  a  more  appropriate  value  of  Rex  could  be  obtained  by  using  the 
integrated  value  along  the  plate.  This  would  tend  to  reduce  the  local 
value  of  Re  compared  to  that  reflected  in  Figure  129,  causing  the  data  to 
rotate  in  a  clockwise  manner;  indicating  somewhat  lower  skin  friction  as 
a  result  of  the  adverse  pressure  gradient. 

*  , 

In  Figure  130  a  similar  comparison  is  shown  for  6  /x  versus  Re  .  As 
can  be  seen  in  this  figure,  both  the  measured  and  theoretical  boundary 
layer  displacement  thickness  growth  rate  decrease  due  to  the  adverse  pressure 

gradient . 

Shown  in  Figure  131  is  the  rate  of  change  of  the  shape  factor  with  Re  . 
The  indicated  shape  factor  decrease  is  primarily  caused  by  the  boundary 
layer  edge  Mach  number  reduction  since  there  is  no  significant  observable 
distortion  in  the  measured  profiles  due  to  the  adverse  pressure  gradient. 

This  is  consistent  with  the  rate  of  growth  of  the  momentum  thickness. 
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The  mass  flow  in  the  boundary  layer 


ni  Versus  x  for  Cowl-Off  Data  - 
BL _ 

is  plotted  in  Figures  132  and  133  versus  the  axial  station.  These  figures 

exhibit  a  significant  increase  in  mass  flow  rate  over  that  expected  for  zero 

pressure  gradient  flow  as  indicated  by  the  theoretical  curves  superimposed  on 

Figure  132.  Figure  132  also  indicates  that  the  trip  configuration  used  for 

the  Model  E-3  test  at  Mach  10,5  did  not  introduce  a  significant  additional 

mass  flux  into  the  boundary  layer  as  evidenced  by  the  agreement  of  the 

tripped  and  untripped  data. 

Cowl-Off  Boundary  Layer  Edge  Mach  Number  Comparison  -  Theoretical  and 
experimental  comparisons  of  the  boundary  layer  edge  Mach  number  distribution 
along  the  ramp  are  shov/n  in  Figure  134  for  the  Mach  10.5  and  7.3  test 
conditions.  The  boundary  layer  edge  Mach  numbers  measured  during  the  E-3 
test  were  somewhat  lower  than  both  the  inviscid  and  viscous  theoretical 
estimates.  Each  data  point  represents  results  form  a  different  run, 
representing  slightly  different  test  conditions  and  thus,  is  also  a  source 
of  data  scatter.  Additional  scatter  is  not  doubt  caused  by  the  fact  that 
both  the  static  pressure  at  the  boundary  layer  edge  and  the  absolute  boundary- 
layer  height  are  difficult  to  establish  more  precisely  from  the  data  due  to 
the  normal  static  pressure  gradient.  For  example,  for  the  low  pressure 
Mach  10.5  data,  an  error  of  .005  psi  in  the  static  pressure  measurement 
represents  an  error  of  .5  in  the  boundary  layer  Mach  number.  (The  repeat¬ 
ability  of  the  pressure  calibration  and  measuring  system  used  during  the 
test  was  equivalent  to  .005  psi.) 

Transition  Study  -  To  aid  in  the  boundary  layer  data  reduction,  the 
power  law  parameters  "n"  (compressible)  and  "n"  (incompressible)  were 
computed  by  the  method  of  least  squares.  The  power  law  parameter  was  again 
used  for  appraising  the  state  of  the  boundary  layer;  laminar,  transitional 
or  turbulent. 

The  results  of  this  analysis  are  also  tabulated  in  Table  III.  As  can 
be  seen  from  this  table,  for  the  high  Reynolds  number  run,  the  boundary  layer 
layer  transitioned  at  approximately  the  same  station  as  the  Model  C  data. 

This  is  not  unusual  in  that  at  these  test  conditions  transition  occurred 
in  a  very  weak  pressure  gradient  field.  For  the  low  pressure  runs,  however, 
transition  occurred  earlier  for  the  Model  E-3  data  compared  to  that  for 
Model  C.  Transition  for  this  case  occurred  in  the  stronger  adverse 
pressure  gradient  field  which  promoted  the  earlier  transition. 

The  transition  Reynolds  number  based  upon  momentum  thickness  was 
calculated  for  several  of  the  transitional  Model  E-3  runs.  Since  boundary 
layer  profiles  were  measured  at  the  three  discrete  probing  stations,  the 
beginning  and  end  of  transition  is  not  precisely  defined.  The  Reynolds 
number  at  these  probe  stations,  however,  was  calculated  and  plotted  on 
Figure  55.  To  indicate  the  degree  of  transition,  the  value  of  the  power 
law  profile  parameter  ("n")  is  also  tabulated  on  this  figure.  Also  shown 
in  the  figure  are  the  data  from  the  Model  A  and  E-l  tests.  From  the  Model  C 
data,  the  transitional  data  is  shown  bracketed  by  the  laminar  and  turbulent 
cases  to  indicate  the  limits  for  boundary  layer  transition.  In  general,  the 


Model  E-3  transition  data  agrees  with  the  Model  E-l  data  when  transition 
occurred  on  the  aft  portion  of  the  model  in  the  strong  adverse  pressure 
gradient  field.  When  transition  occurred  in  the  mild  pressure  gradient  field 
near  the  leading  edge  of  the  model,  the  data  more  closely  agreed  with  the 
Model  A  and  C  flat  plate  data.  The  Model  E-3  trip  data  Indicates  that  the 
trips  used  with  this  model  were  less  effective  in  promoting  earlier  transi¬ 
tion  than  those  trips  generally  employed  on  Model  C. 

In  general,  no  significant  distortion  of  the  profiles  was  observed  for 
e  Model  E-3  data  due  to  the  continuous  adverse  pressure  gradient.  A  slight 
n^stortion  vs 3  indicated  earlier  in  Figure  58  for  the  laminar  boundary  layer 
when  compar  3d  to  the  Blasius  incompressible  profile.  In  Figure  13:?  the 
Model  E-3  turbulent  boundary  layer  is  compared  against  Model  C  profiles  at 
comparable  conditions.  As  can  be  seen  from  this  figure,  no  significant 
distortion  can  be  noted.  Several  profiles  in  the  shock  induced  separation 
region  for  the  extremely  low  Reynolds  number  conditions  were  more  distorted 
as  would  be  expected,  but  much  more  difficult  to  interpret  because  of  severe 
gradients  in  measured  flow  properties. 

Boundary  Layer  Edge  Determination  From  Schlieren  Pictures  -  As  a 
check  on  the  interpretation  of  schlieren  pictures  of  the  boundary  layer, 
composites  of  the  boundary  layer  data  superimposed  on  schlieren  pictures  were 
constructed.  These  composites  are  shown  in  Figures  136  and  137  respectively 
for  typical  laminar  and  turbulent  boundary  layer  cases.  The  outer  boundary 
of  the  apparent  boundary  layer  edge  in  the  schlieren  pictures  for  both  the 
laminar  and  turbulent  cases  corresponds  to  the  real  boundary  layer  edge  based 
on  the  profile  measurements.  The  differences  between  the  laminar  and 
turbulent  boundary  layer  become  very  apparent  when  the  density  profile  is 
plotted  on  the  schlieren.  The  laminar  profile  is  characterized  by  a  sharp 
density  gradient  near  the  outer  edge  of  the  profile  which  corresponds  to 
the  light  line  in  the  schlieren  picture.  The  outer  edge  of  this  line 
defines  the  boundary  layer  edge.  The  turbulent  layer  has  a  more  gradual 
density  gradient  from  the  wall  to  the  outer  edge  of  the  boundary  layer 
thickness.  The  relationship  of  6*  to  the  outer  edge  of  the  boundary  layer 
is  also  shown  in  the  figures.  It  should  be  noted  that  the  light  areas  in 
the  schlieren  could  be  changed  by  adjusting  the  position  of  the  knife  edge 
used  in  the  schlieren  system,  but  the  same  observations  could  be  made,  the 
only  difference  being  the  light  areas  would  be  dark  and  vice  versa.  This 
procedure  for  determining  the  boundary  layer  edge  was  used  in  establishing 
the  boundary  layer  thickness  ratio  across  the  interaction  region  for  Mo^el  C 
in  the  absence  of  desired  profile  traverse  data. 

(2)  Cowl-On  Configuration 

Comparison  of  Theory/Data  Maximum  Pressure  Rise  -  The  maximum  theoretical 
pressure  rise  for  Model  E-3  cowl-on  testing  was  computed  (neglecting  viscous 
effects)  as  a  function  of  cowl  angle  for  pre-test  purposes.  The  procedure 
for  calculating  this  pressure  rise  was  to  assume  that  the  turning  angle 
through  the  shock  from  the  cowl  remained  constant  through  the  compression 
field  generated  by  the  Model  E-3  surface.  The  pressure  rise  across  the 
shock  reflection  was  then  calculated  using  local  flow  conditions  at  the 
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impingement  point  on  the  cuived  ramp.  Comparisons  of  the  data  and  these 
predictions  are  shown  in  Figure  138  and  139  for  the  Mach  7.3  and  10.5  test 
conditions  respectively.  The  data  compares  very  well  with  these  estimates 
for  the  low  cowl  angles.  At  the  high  cowl  angles,  however,  the  data  is 
significantly  lower  than  the  predicted  levels.  This  is  caused  by  the 
expansion  off  the  trailing  edge  of  the  cowl  which  keeps  the  pressure  rise 
from  achieving  its  full  value.  (The  Model  E-3  cowl  was  not  originally 
designed  for  the  high  cowl  angles  subsequently  used  during  this  program. 
During  testing,  thr  cowl  was  set  at  these  high  angles  to  achieve  otherwise 
unobtainable  increased  shock  strengths.) 

Composite  Comparison  of  Pressure  Distribution  as  a  Function  of  Cowl 
Angle  -  A  set  of  composite  pressure  distributions  is  shown  for  the  high  and 
low  Reynolds  number  runs  for  Mach  7.3  and  10.5  in  Figures  140  to  143.  A  set 
of  schlieren  pictures  (Figures  144  to  147)  is  also  shown  for  typical  inter¬ 
action  regions  for  both  laminar  and  turbulent  boundary  layer  at  M  =  7.3  and 
10.5.  As  can  be  seen  from  this  set  of  figures,  no  recognizable  separation 
plateaus  were  observed  at  Mach  10.5  and  7.3  for  test  conditions  corresponding 
to  those  at  which  separation  was  observed  with  Model  C.  Separation  was, 
however,  observed  at  an  extra  low  Reynolds  number  test  condition  at  Mach 
7.3  (Re0  ■  2xl05Ft-^).  This  boundary  layer  stability  trend  is  believed  due 
to  the  fact  tnat  even  at  moderately  low  Reynolds  number  runs,  the  boundary 
layer  transitioned  on  Model  E-3  before  the  shock  interaction  as  a  result 
of  the  extreme  aft  position  of  the  impingement  and  the  fact  that  the 
upstream  adverse  pressure  gradient  promoted  earlier  transition.  Consequently, 
subsequent  testing  was  conducted  at  extremely  low  Reynolds  number  test  condi¬ 
tions  as  stated  above. 

Examination  of  the  pressure  data  for  Runs  65  and  88  (Figure  148) 
indicates  significant  differences  in  spite  of  comparable  shock  strength, 
cowl  position  with  only  a  30%  difference  in  Reynolds  number.  Correspondingly, 
studies  of  the  schlieren  photograph  for  Run  88  (Figure  149)  conform  the 
extensive  separation  region  evident  in  the  pressure  distribution,  whereas 
the  schlieren  picture  for  Run  65  also  shown  in  Figure  149  exhibits  what 
appears  to  be  a  precompression  shock  ahead  of  the  interaction  region.  The 
pressure  distribution,  however,  does  not  indicate  a  typical  pressure  plateau. 
The  mechanism  responsible  for  this  apparent  discrepancy  between  the  two 
otherwise  reasonably  comparable  runs  is  not  evident  from  the  data  obtained. 

Separation  Data  -  As  previously  mentioned  no  laminar  separation  was 
experienced  on  Model  E-3  at  test  conditions  comparable  to  the  Model  C  test 
conditions.  However,  at  Mach  7.3  an  extra  low  Reynolds  number  test  was 
conducted  to  achieve  laminar  separation.  Very  distinct  and  extensive 
separation  plateaus  were  observed  at  these  conditions.  The  plateau 
pressure  levels  for  these  conditions  are  compared  to  predictions  based  on 
the  method  of  Erdos  and  Pallone  in  Figures  150  and  151. 

The  spread  shown  in  the  Erdos  and  Pallone  prediction  is  due  to  the 
uncertainty  in  the  boundary  layer  edge  Mach  number  which  is  due  to  a  possible 
static  pressure  error  of  .005  psi.  In  general,  the  prediction  of  Erdos  and 
Pallone  is  slightly  lower  than  the  measured  pressure  plateau  pressure  for 
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the  separation  data  taken  during  the  Model  E-3  test. 

At  the  highest  cowl  angle  (18°),  turbulent  boundary  layer  separation  was 
detected  during  Mach  7.3  as  indicated  in  the  schlieren  and  pressure  data  of 
Figures  152  and  153.  Although  the  pressure  data  alone  might  be  inadequate 
for  positively  confirming  the  existence  of  a  separation  plateau,  the  pre- 
compression  generated  by  the  separated  flow  region  is  quite  evident  in  the 
schlieren  photograph. 

Summary  -  The  mild  adverse  pressure  gradient  field  had  only  a  minor 
effect  on  the  boundary  layer  properties  as  borne  out  by  direct  comparison 
with  similar  Model  C  flat  plate  profile  data.  In  general,  a  slight  distor¬ 
tion  of  the  laminar  profile  was  noted,  however,  there  was  no  appreciable 
distortion  of  the  turbulent  boundary  layer  profile.  There  was  also  a  slight 
decrease  in  the  measured  laminar  skin  friction  compared  to  Model  C  as  deter¬ 
mined  from  momentum  thickness  data,  while  the  displacement  thickness  and 
shape  factor  decrease  was  consistent  with  the  Mach  number  distribution  for 
the  model.  The  boundary  layer  transition  results  are  consistent  with  the 
Models  A,  E-l  and  C  data.  The  flow  field,  pressure  distribution  and  over¬ 
all  shock  pressure  ratio  were  consistent  with  theory.  No  laminar  boundary 
layer  separation  was  observed  on  Model  E-3  at  conditions  comparable  to 
Model  C.  This  is  believed  due  to  the  fact  that  the  boundary  layer  transi¬ 
tioned  before  the  shock  interaction  as  a  result  of  the  extreme  aft  position 
of  the  shock  impingement  and  the  fact  that  the  upstream  adverse  pressure 
gradient  promoted  the  earlier  transition.  Consequently,  test  Reynolds 
numbers  were  further  reduced  and  laminar  separation  was  then  achieved.  The 
measured  plateau  levels  were  slightly  higher  than  predicted  based  on  condi¬ 
tions  existing  Just  upstream  of  the  interaction.  Turbulent  boundary  layer 
separation  was  noted  at  the  higher  shock  strengths. 

E.  CONCLUSIONS 

As  a  result  of  the  final  phase  of  this  hypersonic  incident  shock-wave 
boundary  layer  interaction  program,  the  following  is  concluded: 

1.  The  modifications  made  to  the  Laminar-Shock/Boundary  Layer  Inter¬ 
action  program  represent  extended  capability  in  obtaining  solutions 
for  both  nonadiabatic  wall  and  adiabatic  wall  laminar  boundary  layer 
incident  shock  problems.  Future  effort  in  this  area  should  be 
directed  at  the  use  of  non-similar  profiles  and  methods  of 
analytically  handling  downstream  transition. 

2.  The  data  from  the  tests  conducted  with  the  3  incident  shock  models 
generally  correlated  well  with  existing  theories  for  boundary  layer 
characteristics  in  zero  and  continuous  adverse  pressure  gradients. 
The  shock  induced  laminar  separation  pressure  plateau  levels  were 
in  excellent  agreement  with  the  theory  of  Erdos  and  Pallone.  The 
laminar  boundary  layer  separation  lengths  were  found  to  increase 
with  increasing  shock  strength,  and  exhibited  an  apparent  Mach 
number  dependence  not  predicted  by  the  Erdos  and  Pallone  formula¬ 
tion.  The  hypersonic  turbulent  boundary  layer  did  not  separate 
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2.  (Continued)  at  conditions  originally  predicted  on  the  basis  of 
lower  Mach  number  theory. 

3.  Analysis  of  the  test  data  aided  in  the  formulation  of  a  flow  model 
for  shock  impingements  on  a  turbulent  boundary  layer.  Solutions 
subsequently  obtained  across  the  immediate  region  of  the  interaction, 
based  on  the  formulated  flow  model,  wore  found  to  be  in  general 
agreement  with  the  test  data  in  .erms  of  boundary  layer  thickness 
change  and  the  prediction  of  imminent  separation.  The  strong  effect 
of  upstream  profile  shape  and  wall  to  total  temperature  ratio  on 

the  separation  characteristics  predicted  by  the  solutions  obtained, 
delineate  the  importance  of  these  two  parameters  with  respect  to 
future  experimental  programs. 

4.  The  technique  developed  for  estimating  the  change  in  turbulent 
boundary  layer  properties  across  the  immediate  region  of  a  shock 
impingement  is  adequate  for  use  in  current  engineering  estimates 
associated  with  the  design  and  analysis  of  hypersonic  inlets.  Sub¬ 
sequent  effort  in  this  area  should  be  directed  at  extending  the 
analysis  to  include  downstream  constant  pressure  mixing,  considering 
the  effects  of  initial  flow  profile  distortion. 

5.  For  these  hypersonic  inlet  applications  where  the  laminar  boundary 
layer  has  not  transitioned  by  the  first  shock  impingement  point, 
shock  induced  transition  should  be  expected  shortly  downstream.  For 
shock  strengths  in  excess  of  4°  turning,  laminar  separation  should 
be  expected.  Heat  transfer  estimates  downstream  of  the  shock  im¬ 
pingement  point  should  reflect  appropriate  turbulent  boundary  layer 
heating  rates. 

6.  Reduced  wall  to  total  temperature  ratios  associated  with  full  scale 
inlet  operation  at  flight  enthalpy  levels  represent  a  potential 
source  of  increased  turbulent  boundary  layer  stability  compared  to 
separation  characteristics  at  higher  wall  to  total  temperature 
ratios.  This  should  be  considered  in  the  interpretation  of  inlet 
component  test  data  commonly  obtained  in  gr <und  test  facilities  at 
higher  wall  to  total  temperature  ratios  when  full  scale  enthalpy 
levels  are  not  simulated.  Assuming  attendant  thermal  protection 
requirements  can  be  met,  flight  enthalpy  hypersonic  inlet  boundary 
layer  separation  problems  may  be  alleviated  through  the  use  of 
appropriately  designed  boundary  layer  trips. 
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SECTION  IV 
APPENDIX  I 


Laminar-Shock/ Boundary  Layer  Interaction  Program  Users' 

Manual 

a.  Program  Description 

The  Laminar-Shock/Boundary  Layer  Interaction  Program  is 
intended  to  calculate  the  wall  pressure  distribution,  the  boundary  layer 
growth,  and  the  wall  shear  in  a  region  of  sudden  pressure  rise  and  pos¬ 
sible  boundary  layer  separation.  The  principal  application  is  that  of  a 
shockwave,  induced  by  a  secondary  surface,  impinging  on  a  flat  plate 
or  wedge  in  uniform  supersonic  flow.  A  second  application  is  the  boun¬ 
dary  layer  inviscid  flow  interaction  in  a  concave  corner,  as  formed  by  a 
double  wedge.  The  concave  corner  (with  an  exterior  angle  of  2A  9)  pro¬ 
duces  an  effect  similar  to  an  impinging  shock  wave  (with  flow  deflection 
A0).  The  program  is  limited  to  small  surface  inclinations,  to  a  laminar 
boundary  layer,  and  to  a  single  right  running  (impinging)  shock  wave. 

Right  characteristics  (exclusive  of  the  shock  wave)  must  be  of  zero 
strength.  The  program  has  some  limited  capability  for  handling  curved 
flow  boundaries. 

The  calculation  procedure  is  based  on  a  single  analytical  model 
in  which  the  separated  and  attached  boundary  layer  flows  are  treated  with¬ 
out  distinction.  In  both  regions  the  boundary  layer  equations  are  integrated 
simultaneously  with  the  Prandtl-Meye r  inviscid  flow  relations.  A  concise 
statement  of  the  problem,  the  formulated  equations,  and  a  brief  description 
of  the  solution  procedure  is  presented  in  this  User's  Manual. 


b.  Outline  of  Shock-Boundary  Layer  Interaction  Calculation 


1.  Surface  coordinates,  Yg 

2.  Shock  wave  incidence  point,  X2  =  Xg^, 
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3. 


Shock  wave  strength,  i.e.  ,  flow  deflection  through  shock 
wave  ^  9 


Starting  conditions  at  Xj 
Mej.  n,,  and  H,r 


1 


Conditions  at  X4 

Me^,  n^,  and  Htr^  (Me^  is  determined  by  Me  ^  and 

the  shock  strength) 


Major  Assumptions: 

1.  Perfect  gas  with  constant  C  . 

6  p 

2.  Viscosity  proportional  to  temperature 

3.  Unity  Prandtl  number 

4.  Two-dimensional  flow 

5.  Prandtl-Meyer  inviscid  flow  up  to  and  downstream 
of  the  shock  incidence  point, 

6.  The  incident  shock  wave  is  reflected  from  the  boundary 
layer  edge  as  from  a  constant  pressure  boundary. 

7.  Zero  Pressure  gradient  across  the  boundary  layer, 


Unknowns : 

The  unknowns  are  divided  into  three  categories  as  follows: 

1.  Xj 

2.  X^  and  84 

3.  9,  Me,  n  and  Htr  between  Xj  and  X4 

(Xj  is  the  extent  of  forward  feed  of  the  shock  induced  pressure  rise, 
X4,  denotes  the  termination  point  of  the  interaction  zone) 

Method  for  Integration  of  Boundary  Layer  Equations 
(Solution  for  unknowns  in  category  3): 
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Trial  values  of  Xj,  X4  and  64  are  assigned.  Then  the  boundary 
layer  equations  are  integrated  between  the  calculation  starting  point, 
Xj  or  X4,  and  the  shock  incidence  point,  X2.  Compatibility  between 
the  boundary  layer  growth  and  the  Prandtl-Meyer  inviscid  flow  is 
enforced. 


Isentropic  flow  is  assumed, so  that  the  boundary  layer  edge  condi¬ 
tions  (namely  the  quantities  aeP  and  and  the  Prandtl-Meyer 

angle  are  functions  of  the  boundary  layer  edge  Mach  number  Me.  The 
boundary  layer  parameters  L ,  E/0  G,  (6-6*)/9,  and  cp  are  functions  of 
Htr.  The  equations  are: 


1. 


Von  Karman  momentum  integral 


_d 

dx 


( aeP) 


I  Pe  ae  9^ 

1 


2 


aeP 

Me 


1  +  n  (Htr  +  2) 


(Al) 


2.  Equation  defining  pressure  gradient  parameter 


dMe 
d  x 


n 


(  Pea.82 

l 


3. 


Kinetic  ene rgy/momentum  integral 

n  =  i^-Gt 

G  ( Htr  -  1)  -  2E 

9 


(A2) 


(A3) 


4. 


Boundary  layer  growth-inviscid  flow  compatibility  (v 
determined  by  initial  conditions  at  X]  or  X4). 


6  =  0 


(Htr 


Y-  1 


+  1)  ( 1  +  2 


Mp2) 


-  1 


ref 


is 


(A4) 


D  = 


Me2  -  1  2 

1  +  Me 

2 


n 


Ppue 


) 


(A5) 
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The  above  equations  are  numerically  integrated  by  employing  the 
trapezoidal  rule  over  a  small  interval  defined  by  ZX  -  Xj  -  Xj_  1 .  The 
integration  proceeds  as  follows: 


1.  Assume  values  of  H*.-.,  m  and  M...  at  X;,  the  end  of  the 

iri  1  ci 

interval. 

2.  Determine  pe  ae02  by  integration  of  Equation  (Al). 

£ 

Notice  that  aeP  is  a  function  of  Me 

3.  Determine  Me-  by  integration  of  Equation  (A2)  . 

4.  Repeat  steps  2  and  3  until  agreement  is  obtained  between  the 
assumed  and  calculated  MP.. 

5.  Calculate  ni  by  Equation  (A3) .  Note  that  dG/dX  is  based  on 
the  assumed  variation  of  Htr  over  the  interval. 

6.  Repeat  steps  2  through  5  until  agreement  is  obtained  between 
the  assumed  and  calculated  nj. 

7.  Evaluate  both  sides  of  equation  (A6) .  If  agreement  is  not 
obtained,  assume  a  new  value  of  Htr,  and  repeat  steps  2 
through  6. 

Iteration  Procedure  for  Xj,  X4  and  64: 

Typical  results  from  the  integration  of  the  BL-PM  equations 
using  trial  values  of  X],  X4  and  84  appear  as  follows: 
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Compare  Me^  with  Me^ 


Agree  Do  not  agree 


Choose  new  X  j 

Final  Solution  Obtained  j 
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c.  Input 


A  table  of  boundary  layer  properties  is  first  read  into  the  program 
under  the  NAMELIST  name  NSSBL;  presently  tables  are  available  for 
two  wall  temperature  ratios.  Any  number  of  sets  of  input  for  specific 
problems  are  then  read  under  the  NAMELIST  name  SH0CKB.  The  input 
can  be  prepared  and  keypunched  from  input  sheets  as  illustrated  by  the 
example  on  page  6l‘.  The  input  quantities  are  defined  on  page  59. 
Additional  comments  regarding  the  input  items  are  as  follows: 

The  Mach  number  and  Reynolds  number  are  input  for  flow  conditions 
parallel  to  the  x-axis.  For  example,  if  the  flow  surface  is  a  3°  wedge 
we  would  let  AA  =  CC-0  and  BB  =  tan  3°.  For  a  first  approximation  the 
free  stream  Mach  number  and  Reynolds  number  could  be  used  for  Q  and 
RE1.  However,  by  reference  to  Figures  154and  155Q  and  RE]  can  be 
easily  corrected  to  account  for  the(average)  pressure  recovery  across 
the  leading  edge  shock  wave.  In  these  Figures  the  Mach  number  and 
Reynolds  Number  ratios  are  plotted  vs.  total  pressure  recovery  assuming 
the  flow  is  re-expanded  to  the  axial  direction. 

The  impinging  shockwave  strength  is  input  by  specifying  the  flow 
deflection  through  the  wave,  D0S.  D0S  is  zero  for  the  concave  corner 
problem. 


The  integration  step  size  is  determined  by  DXX.  For  most  con¬ 
ditions  the  following  relation  may  be  used  to  estimate  an  appropriate 
value  of  DXX. 


DXX  -  .  25 


(2.  59  TTe  + 
TW 


l)  ( 1  +  m2) 

2 


1 


_ X _ 

.  1/4  1/2 

pux  I  (B) 


BN,  B1  and  B4  are  pressure  gradient  parameters.  BN  is  only  used 
in  the  calculation  of  momentum  thickness  at  the  upstream  boundary  XI. 
Normally  this  value  would  be  zero;  it  may  be  input  to  duplicate  the  value 
of  BN  obtained  from  a  leading  edge  interaction  calculation. 


BN 


X_  dP 
P  dx 


B1  and  B4  are  actual  initial  conditions  of  the  parameter  n  (for  the 
boundary  layer  integrations)  at  the  upstream  and  downstream  stations 
1  and  4  where  n  is  defined  as: 
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n 


feae  6  dMe 

Me  dx 

The  relation  between  BN  and  n  is  shown  in  Figurel56  .  However,  in 
order  to  obtain  a  meaningful  interaction  solution,  the  value  of  B1  must 
be  arbitrarily  set  to  a  value  greater  (more  positive)  than  the  precise 
value  obtained  from  Figurei56  for  known  values  of  BN.  Thus,  B1  is 
usually  set  to  zero,  although  BN  and  the  corresponding  n  are  slightly 
negative . 

Similar  rules  must  be  followed  at  the  downstream  boundary, 
station  4.  Here  a  value  B4  =  .01  is  often  used  to  artificially  reduce 
the  distance  to  the  terminal  point  X4;  the  results  in  the  immediate 
vicinity  of  the  shock  incidence  are  only  very  slightly  affected. 

A  sample  input  sheet,  prepared  for  the  problem  of  Figure  3,  is 
presented  on  page  61.  The  input  sheet  itself  is  that  currently  used 
with  the  computer  program. 


d.  Sample  Output 

A  sample  output  is  shown  in  Figure  157  (see  corresponding  plot 
in  Figure  3).  Parameters  relating  to  the  boundary  layer  growth  and 
external  pressure  distribution  are  printed.  A  list  of  the  output  symbols 
follows : 

Input 


AA,  BB,  &  CC 

Coefficients  for  quadratic  surface.  See 
input  sheet. 

BN,  B 1  &  B4 

Defined  in  Section  3. 

FSN 

Weighting  factor;  see  input  sheet. 

X  (SHOCK) 

Axial  position  of  shock  wave  impingement 

or  convex  corne  r. 

H(W)/H(TOT) 

Input  wall  to  stagnation  enthalpy  ratio. 

DEL  THETA 

Flow  deflection  through  shock  wave,  deg. 

REY 

Effective  unit  Reynolds  number  corresponding 
to  Mach  number,  MACH. 

MACH 

Mach  number  of  flow  in  axial  direction. 
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GAMMA 


Ratio  of  specific  heats 


Output 

DDSDX 

X 

N  (HAT) 


MACH 

P/Pl 

N 

\ 


HTR 

L 

PHI 

THETA 
DELTA* 
INV.  Y 
WBL 


Station  1  value  of  d  6* 

dx 


Axial  coordinate 

Second  derivation  parameter 


n  -  _9_^_  _H£_  dT  1 

Meue  he  [  d  Y  J  w 

Boundary  layer  edge  Mach  number 

Pressure  ratio;  PI  pressure  corresponding 
to  the  input  MACH  number. 

Pressure  gradient  parameter 

n  =  -  Pe6^  He  due 

14;  he  dx 


e2  hq  dp 

^ue  he  dx 


Low  speed  form  factor,  Ht 
Shear  parameter, 

Kinetic  ene  rgy  dis  sipation  function, 
Momentum  thickness,  0 

* 

Displacement  thickness  ,  6 
Ordinate  of  boundary  layer  edge. 
Boundary  layer  flow  parameter 


(  6  -  6*  ) 


0e«e 


(P0U0)  « 


Subscript  "o"  refers  to  conditions  at  the 
input  MACH  number. 
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Page _ of 


SHOCK-BOUNDARY  LAYER  INTERACTION 
F0988A  SHOCK-BL  NAMELIST 


Plus  signs  may  be  omitted.  Use  exponents  for  long  numbers. 
Enter  30  X  10°  as  30E6.  Enter  .  15  X  10^  as  .  15E6.  Do  not  use 
commas  within  field,  i.  e.  ,  enter  3, 800  as  3800. 


required  input 

W 

$ 

$  SHOCKB 

ident  -6oh  0  sTeaMr-.fr>  sHocffo  cohJvErr,EklC£ 
Q-  10.4  ,  RE  1-  L.32&4  ,  XSK=  M  D8S=  5  ,  GAMMA- 

surface  specifications 
option  A 

DXX-  |  ,  X  =  \Z  AA^  O  BB=.q524I  ,  CC^  Q  , 

option  B 


f 

9 

1 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

i 

1 

f 

9 

9 

9 

9 

f 

_ 9  _ 

_ 9  _ 

9 

_ _ 9  _ 

_ 9  _ 

_ _ 9 

K= _ ,  NK  - 


optional  input 

X4I= _ ,  941= _ ,  D03D94  = 

BN=  -  2  ,  BI=  Q  ,  B4=  ,  o  / 

FNS= _ , 

STOP= 
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NOTE:  Input  also  includes  a  table  of  boundary  layer  properties  for  the 

appropriate  wall  to  stagnation  temperature  ratio.  Tables  are 
presently  available  for  T^/Tj  0.  2  and  1.  0. 


SHOCK  -BL 

EXPLANATION  OF  TERMS 


1 .  Required  Input: 

Q  Mach  number  of  flow  achieved  by  isentropic  expansion  (or 
compression)  to  the  axial  direction. 

RE  1  Effective  unit  Reynolds  number  corresponding  to  Mach  number 
Q  (including  non-linear  viscosity  correction),  RE1  p  u 

C 

^o  Tw 

XSK  Axial  position  of  shock  wave  impingement  or  convex  corner 


D0S  Flow  deflection  through  the  shock  wave,  deg. 

GAMMA  Presently,  the  deck  is  limited  to  either  y  1.4  or 

V  =  1.667.  If  not  specified,  y  1.4  will  be  assumed, 

2 .  Surface  Specification: 

Option  A 

DXX  Axial  interval  of  the  integration  steps  times  four 

X  First  guess  of  X  1 ,  the  beginning  of  the  interaction  region 

AA,  BB,  CC,  Coefficients  for  quadractit  surface* 

y B  AA  4  BR'X  4  CC-X2 

to  be  omitted  or  set  equal  to  zero  for  nominal 
horizontal  surface. 
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Option  B 

X,  Y  List  of  surface  coordinates  -  maximum  of  100.  Note  that  the 
integration  interval  is  (Xn  -  Xn_j)  /4.  Note  that  for  an  arbitrary 
curved  boundary,  the  coordinates  must  be  specified  close  together 
with  great  accuracy. 

K  First  guess  of  XI  (beginning  of  interaction  region)  is  Kth  value 
of  X. 

NK  Number  of  X-Y  pairs  input 

3.  Optional  Input  for  Iteration 

X41  First  guess  of  axial  position  of  station  4,  the  terminal  point 
of  the  interaction  zone. 

041  First  guess  of  momentum  thickness  at  station  4. 

4.  Optional  Input  -  Unless  other  wise  specified,  zero  is  used  for  the 

following  parameters. 

BN  Power  law  exponent  for  pressure  variation  upstream  of  XI 
fp  a  j  used  in  the  calculation  of  momentum  thickness 

at  XI.  The  boundary  layer  is  assumed  to  start  at  X  -  0. 

Bl,  B4  Boundary  conditions  on  n  at  stations  1  and  4. 

FSN  Weighting  factor  on  the  use  of  wall  momentum  condition 
(FSN  =  1)  or  K.E.  equation  (FSN  =  0) 

STOP  Specify  STOP  =  2  for  trial  solution  without  iteration  for 

XI,  X4  &  04.  Specify  STOP  =  2  for  upstream  free  interaction 
calculation  only. 

D03D04  Rate  of  change  of  8  at  3  compared  to  rate  of  change  of  9  at  4. 
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APPENDIX  II 


Using  the  computational  procedure  for  estimating  the  change  in 
turbulent  boundary  layer  thickness  across  an  impinging  shock  wave  using 
the  equations  formulated  in  Reference  5,  the  boundary  layer  thickness 
change  across  a  turbulent  boundary  layer-shock  wave  interaction. 


can  be  estimated  as  follows: 

Tw  ,  P3 

1.  Given  Mj,  M3,  -^7—  -p— 

(compressible)  ^ 


and  upstream  profile  properties 


2.  Compute  upstream  incompressible  profile  parameter: 
Nl  =  Z(  TW/TTe)  (1+  ) _ 


H  -  Me2-(TW/TTe)  (1  +  Me2) 


3.  (Hj),  = 


Ni  +  2 


N 


1 


4. 


f(Hi) 


1.0U 


-IdHi 

- or  from 


Hj  ( Hj2  -  1)  (Hj  +  1) 

1- 

[l-TW  1 

(Hi2+4Hi-l) 

1  tt  1 

1  e 

(Hi+l)(Hi+3) 

m 

curves  of  Figure  158. 
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5.  flH;),  =  In  ^1+  «H;)1 

Mj 

6 .  Solve  (4)  for  ( 3 - >  if  (Hi)3  >3  — ►  SEPARATION 

7.  N3  =  2 _ 

"  ( Hi ) 3  -  1 

5 6* 

8.  At  Mj,  Ni  and  M3>  N3  compute  (—  )j  and  (~)3 


Tw  (  1  N 

Txe  N+l  1+  M  2  (N+l)(N+2) 
_  2  e  _ 


TW  +  ( i  .  TWj  N 
T 'j’  T N+ 1 


N 

N+2 


1  -  (6*/6)3  M3 

1  .  (  6=:=/  6)j  Mj 


10.  6*3  _  (6*/6)3  l 

~6^  (  6*8*/ 6)j  (  6  j  /  6  3 ) 
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appendix  III 


Data  Reduction  Techniques 
a.  Boundary  Layer  Data  Reduction 

Introduction 

To  expedite  the  analysis  of  the  large  amount  of  profile  data 
obtained  during  the  Model  C,  D  and  E-3  testing  program,  a  computer 
program  was  utilized  which  would  take  as  input  the  measured  boundary 
layer  data  as  tabulated  in  Reference  11  and  Reference  12,  which  are  the 
Model  C,  D  and  E-3  data  reports. 

The  computer  data  reduction  program  can  be  run  with  any  of 
four  options  on  the  selection  of  input  data.  The  difference  between  these 
four  options  is  as  follows: 

1.  Option  I  -  All  the  measured  boundary  layer  data:  probe 
location,  pitot  pressure,  static  pressure  and  total 
temperature  are  used  in  evaluating  the  boundary  layer 
flow  properties. 

2.  Option  II  -  Use  is  made  of  an  assumed  constant  static 
pressure  across  the  boundary  layer  rather  than 
measured  values. 

3.  Option  III  -  A  Crocco  temperature  profile  is  used  instead 
of  the  measured  temperature  distribution,  however,  either 
the  measured  static  pressure  or  a  constant  static  pressure 
can  be  used.  In  either  case,  the  boundary  layer  edge 
velocity  is  assumed  equal  to  freestream  value  in  the 
iteration  required  for  determining  the  velocity  distribution 
within  the  boundary  layer. 

4.  Option  IV  -  This  option  is  similar  to  Option  III  with  the 
exception  that  the  boundary  layer  edge  velocity  as 
determined  from  a  previous  pass  using  Option  III  is  used 
in  the  total  temperature  calculations  instead  of  the  free¬ 
stream  velocity. 
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The  equations,  assumptions  and  procedures  used  in  calculating 
the  boundary  layer  properties  with  these  different  options  are  described 
in  the  next  subsection. 

For  nearly  all  of  the  wedge  Model  C-D  data  reduction,  Option  IV 
was  used  with  constant  static  pressure  assumed  across  the  boundary  layer. 
For  the  curved  surface  Model  E-3  data  reduction,  Option  IV  was  used  with 
variable  static  pressure  across  the  boundary  layer. 

b.  Data  Reduction  Procedure 

Boundary  layer  surveys  were  made  using  the  NASA  boundary 
layer  probe  of  Figure  31.  The  probe  and  its  drive  unit  are  described  in 
Section  III-A  which  also  identifies  the  location  of  the  probing  stations  on 
each  of  the  three  test  models.  Pitot  pressure,  static  pressure,  and 
total  temperature  measurements  were  made  across  the  boundary  layer. 

The  flow  properties  in  the  boundary  layer  were  then  calculated 
on  the  IBM  7094  computer  using  the  selected  Option  and  the  probe 
measurements  as  described  in  the  following  paragraphs. 

The  Mach  number  at  a  point  in  the  boundary  layer  was  computed 
from  the  pitot  and  static  pressure  using  the  Rayleigh  pitot  formula  for 
perfect  gas.  Initially  this  equation  was  corrected  for  caloric  imperfec¬ 
tions  using  the  correction  factor  given  in  Chart  17  (supersonic)  and 
Chart  1 Z  (  subsonic)  of  NACA  Report  1135.  Since  this  corresponded  to  a 
one  percent  or  less  correction  to  the  Mach  number  at  the  Ames  testing 
conditions;  this  correction  was  subsequently  neglected. 

From  the  Mach  number  and  the  total  temperature,  the  static 
temperature  was  determined  using  the  following  perfect  gas  relation. 

,  7-1  ..2 

T  =  Tx/(  1  +  —  M  ) 

Whe  re : 

T  =  static  temperature 
Tp  =  total  temperature 
M  =  Mach  number 
V  -  ratio  specific  heats  =  1.4 
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This  equation  was  corrected  by  the  correction  factor  given  in  Chart  10  of 
NACA  Report  1135.  The  procedure  for  using  the  NACA  1135  correction 
factors  in  the  computer  program  was  to  carpet-plot  the  charts  from  NACA 
1135  and  then  construct  two-dimensional  tables  of  the  correction  factors 
versus  Mach  number  and  total  temperature.  A  two-dimensional  table-look 
up  method,  which  initially  used  parabolic  interpolation,  was  used  to 
determine  the  correction  factor  given  the  Mach  number  and  total  temper¬ 
ature.  The  parabolic  interpolation  routine  was  later  replaced  with  a 
linear  interpolation  routine  to  save  on  computer  time. 

From  the  static  temperature  and  Mach  number,  the  velocity 
was  computed  using  the  definitions  of  Mach  number  and  speed  of  sound 
as  indicated  by  the  following  relation: 

u  =  V  yRT  M 

Where: 


u  = 

velocity,  ft/sec 

R  = 

universal  gas  constant 

T  = 

static  temperature,  °R 

M  = 

Mach  number 

V  — 

ratio  of  specific  heats  at  T 

At  this  static  temperature  the  ratio  of  specific  heats  (  V)  was 
calculated  from  a  curve  fit  of  temperature  versus  Y  as  determined  from 
Keenan  and  Kaye  Gas  Tables. 

The  density  was  then  calculated  from  the  equation  of  state 
using  the  static  pressure  and  temperature.  The  unit  mass  flow  could 
then  be  determined  from  the  product  of  density  and  velocity. 

Using  these  fluid  properties,  the  boundary  layer  integral 
properties  were  then  calculated  by  numerical  integration  across  the 
boundary  layer  using  the  trapezoidal  rule.  The  displacement,  momentum 
and  energy  thicknesses  were  thus  computed  using  numerical  integration 
and  the  following  relations: 


oPeue 


(1  -  T7-  )  dy 


u£ 


JO  p  U 

e  e 


c  u  u 

( i  — t  )  dy 


u 


Where : 
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e 


6 


0 

c 

u 

u 


e 

e 


displacement  thickness 
momentum  thickness 
energy  thickness 
local  density 

boundary  layer  edge  density 
local  velocity 

boundary  layer  edge  velocity 


The  boundary  layer  mass  flow  and  transformed  y  coordinate 
across  the  boundary  layer  are  also  computed.  The  transformed  y 
coordinate  (  p)  is  evaluated  using  the  Stewartson's  transformation: 


■n 


dy 


r)  is  subsequently  used  to  determine  the  incompressible  form  factor  for 
the  boundary  layer. 

For  the  turbulent  boundary  layer  the  incompressible  form 
factor  was  determined  by  plotting  the  velocity  versus  r]  and  measuring 
the  value  of  the  slope.  This  slope  corresponds  to  N  (the  power-law 
profile  parameter)  which  is  related  to  the  incompressible  shape  factor 
by  the  following  relation: 

2 

Hi  -  1  +  F 

Where: 

Hi 
N 


incompressible  shape  factor 
incompressible  power  law  profile  parameter 
in  the  expression  ^ 
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For  the  laminar  cases,  the  r\  transformation  allows  the  data  to  be  compared 
directly  to  the  various  theoretical  results.  (For  example,  see  Figure  159 
for  a  typical  comparison  of  the  experimental  data  to  the  theory  of  Cohen 
and  Reshotko) . 

All  of  the  computed  flow  properties  are  tabulated  at  each  input 
point  along  with  the  integral  properties,  mass  flow  and  r).  Also  tabulated  are 
the  integrands  of  the  various  integral  properties;  the  freestream  and 
boundary  layer  edge  flow  conditions.  Several  options  are  also  available 
for  printing  out  the  data.  One  option  prints  the  boundary  layer  flow 
properties,  another  the  flow  properties  and  all  the  integral  properties, 
while  a  third  prints  the  flow  properties,  selected  integral  properties, 
mass  flow  and  r\.  A  typical  print  out  is  presented  in  Figure  160. 

The  normal  calculational  procedure  was  to  make  a  pass  through 
the  data  using  Option  III  (See  Section  IV-B)  to  determine  the  boundary  layer 
edge  conditions.  Re-inputing  this  data  into  the  computer  program 
Option  IV  was  then  used  to  calculate  the  final  boundary  layer  properties. 

The  reason  for  using  a  two-pass  procedure  was  that  a  check  could  be  made 
on  the  boundary  layer  edge  conditions  before  the  final  integration  was 
conducted. 

c.  Model  C  Data  Reduction 

For  the  Model  C  and  D  boundary  layer  data  reduction,  the 
static  pressure  distribution  was  assumed  to  be  constant  across  the 
boundary  layer  and  equal  to  the  wall  static  pressure  at  the  probing 
station.  The  validity  of  this  assumption  is  verified  by  essentially  constant 
static  pressure  distributions  measured  across  the  boundary  layer  during 
the  Model  C  and  D  testing  program,  such  as  the  typical  example  of 
Figure  161. 

The  total  temperature  measurements  for  Model  C  appeared  to 
never  reach  equilibrium  because  of  a  time  lag  problem  and,  therefore,  a 
Crocco's  total  temperature  distribution  was  used  in  the  data  reduction 
instead  of  the  measured  values. 

d.  Model  E-3  Data  Reduction 

A  re-examination  of  the  boundary  layer  properties  was 
conducted  to  determine  if  the  following  assumptions  were  valid  for 
Model  E-3.  The  assumptions  are: 
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1.  Constant  static  pressure  across  the  boundary  layer. 


2.  A  Crocco  total  temperature  distribution. 

For  the  Model  E-3  data,  the  static  pressure  did  vary  across 
the  boundary  layer  such  that  the  constant  static  pressure  assumption 
could  not  be  used.  A  typical  variation  of  static  pressure  across  the 
boundary  layer  at  the  three  probing  stations  is  shown  in  Figures  162 
through  i64.  As  indicated  in  the  figures,  the  static  pressure  variation 
becomes  more  severe  at  the  aft  probing  position  as  would  be  expected 
since  these  stations  are  in  a  more  severe  axial  adverse  pressure  gradient 
field. 


For  the  low  pressure  runs  {Figure  165),  the  static  pressure 
probe  seemed  to  consistently  give  high  readings  when  compared  to  the 
wall  pressure  measurements.  This  was  possibly  caused  by  increased 
boundary  layer  growth  on  the  probe  at  the  low  Reynolds  number  test 
conditions.  As  a  consequence,  the  static  pressure  distribution  was 
determined  by  matching  the  pressure  at  the  wall  and  using  the  pressure 
variation  measured  by  the  static  pressure  probe. 

Since  it  is  necessary  to  iterate  between  the  velocity  distribution 
and  the  assumed  temperature  distribution,  a  distribution  where  the 
temperature  and  velocity  are  coupled  is  a  desirable  and  necessary  condi¬ 
tion.  Not  enough  information  is  ordinarily  obtained  during  a  boundary 
layer  survey  to  solve  the  momentum  and  energy  equations  at  the  probing 
station. 


The  Crocco  distribution  is  expressed  as 
H  =  Hw  -  (He  -  u/ue 

Where : 

H  =  local  total  enthalpy-BTU/#  mass 
=  wall  enthalpy-BTU/#  mass 
He  =  boundary  layer  edge  total  enthalpy-BTU/#mass 
u  =  local  velocity- ft/ sec 
ue  =  boundary  layer  edge  velocity- ft/ sec 
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The  modified  form  of  the  Crocco  distribution  is  given  as: 

H  =  Hw  +  (He  -  Hw)  +  (R-l)  (u  ue  -  u  )  /2gJ 

Whe  re: 

R  =  recovery  factor 
R  %  \/Pr  laminar,  ^v^Pr  turbulent 
g  =  gravitational  constant  -32.  174  ft/sec 
J  =  Joule's  Constant  -778  ft-#/BTU 

When  Pr  <  1,  the  modified  Crocco  distribution  approximates  the 
enthalpy  distribution  near  the  wall  better  than  the  Crocco  distribution. 

The  modified  Crocco  distribution,  however,  does  not  satisfy  all  conditions 
at  the  boundary  layer  edge;  namely,  a  thicker  thermal  boundary  layer  and 
temperature  overheat  phenomena  occur. 

As  can  be  seen,  the  modified  form  reduces  to  the  standard  form 
when  R  =  1.  A  comparison  of  the  two  theoretical  distributions  are  shown 
compared  with  experimental  data  from  a  typical  high  and  low  pressure  run 
in  Figures  166  and  167  (Run  numbers  10  and  13  respectively  from  the  E-3 
data).  Also  for  Runs  10  and  13  a  comparison  was  made  to  determine  the 
effect,  on  the  integral  properties,  of  using  these  two  temperature  distri¬ 
butions.  There  was  less  than  a  one  percent  effect  on  displacement  thick¬ 
ness  and  approximately  a  five  to  ten  percent  effect  on  momentum  thickness, 
6.  The  overall  effect  was  the  Crocco  distribution  predicted  a  slightly 
fulle r  profile . 

In  general,  the  modified  Crocco  distribution  agreed  better  with 
experimental  temperature  data  near  the  wall  while  the  Crocco  distribu¬ 
tion  gave  better  agreement  near  the  boundary  layer  edge. 

It  was  evident  that  the  temperature  time-lag  problem  still  existed 
during  the  E-3  test  and  can  be  illustrated  in  Figure  168  during  the  up  and 
down  traverse  of  the  total  temperature  prcbe.  In  general,  the  total 
temperature  probe  also  measured  a  total  temperature  at  the  boundary 
layer  edge  approximately  200-300 °R  lower  than  the  tunnel  total  tempera¬ 
ture  (See  Figure  168).  Because  of  this  apparent  time-lag  problem  coupled 
with  the  low  readings,  the  experimental  temperature  data  was  not  used  in 
the  E-3  data  reduction,  but  rather  the  theoretical  distributions  were  used. 
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During  the  Model  C  data  reduction  the  Prandtl  number  was 
assumed  equal  to  one  and  the  Crocco  distribution  was  used  exclusively. 
For  Model  E-3  the  data  reduction  program  was  modified  so  that  either 
temperature  distribution  could  be  used  by  simply  specifying  an 
appropriate  temperature  recovery  factor.  However,  to  be  consistent 
with  the  Model  C  data  reduction,  a  recovery  factor  of  one  was  used  for 
most  of  the  E-3  data. 

e.  Pressure  Distribution  Analysis 

The  data  reduction  for  the  Model  C,  D,  and  E-3  surface  static 
pressure  was  done  at  Ames  by  NASA.  A  considerable  amount  of  data 
was  processed  by  the  NASA  personnel,  nearly  one-half  million  lines. 

In  order  to  reduce  this  data  to  a  less  voluminous  and  more  convenient 
form,  a  condensed  data  report  was  prepared  for  the  Model  C  and  D 
data  and  another  for  the  Model  E-3  (See  References  11  and  12).  The 
data  tabulated  in  these  reports  consists  of  the  raw  data  which  has  not 
been  corrected  for  any  reference  pressure  error,  instrumentation 
time-lag,  etc. 

A  thorough  examination  of  the  data  indicated  that  a  reference 
pressure  error  existed  in  the  data.  It  appears  that  this  error  was 
caused  by  the  calibration  procedure  for  the  pressure  transducers.  This 
procedure  was  to  pump  the  test  cabin  pressure  down  to  .  2  psia  before 
the  test.  The  pressure  transducers  were  then  zeroed  at  this  pressure 
level.  An  error  of  only  2  1/2%  (.005  psia)  in  the  reference  pressure, 
which  is  the  repeatability  of  the  reference  pressure  measuring  system, 
would  cause  approximately  a  20%  error  in  the  surface  pressure  measure¬ 
ments  at  the  low  Reynolds  number,  Mach  10.  5  testing  conditions. 

A  comparison  of  the  surface  pressure  data  and  boundary 
layer  data  against  the  theoretical  expected  values  indicated  that  the 
problem  was  of  a  zero  shift  nature.  In  the  absence  of  any  better  way  for 
predicting  this  reference  pressure  error,  the  pressure  distributions 
were  shifted  to  match  the  G.E./ Bertram  prediction  which  from  past 
experience,  at  the  same  testing  conditions  (for  example  earlier  results 
obtained  under  this  program),  gives  an  excellent  surface  pressure 
prediction.  Therefore,  the  procedure  used  in  the  data  reduction  for 
the  Model  C,  D  and  E-3  test  results  was  to  adjust  the  static  pressure 
readings  by  an  amount  equal  to  the  difference  between  the  leading  edge 
pressure  measurements  and  the  G.E. /Bertram  pressure  predictions. 
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Figure  79  SCHLIEREN-PRESSURE  DISTRIBUTION  COMPOSITE 
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Figure  99  THEORETICAL  SHEAR  PARAMETER 
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igure  143  TYPICAL  TURBULENT  FLOW  AXIAL  PRESSURE  DISTRIBUTION 
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igure  160  TYPICAL  DATA  REDUCTION  PRINTOUT 


P  -  Static  Pressure  -  psia 

Figure  161  TYPICAL  STATIC  PRESSURE  PROFILE 


nches 


F  -  Static  Pressure  -  psia 

Figure  162  TYPICAL  STATIC  PRESSURE  PROFILE 


Inches 


Figure  163 


TYPICAL  STATIC  PRESSURE  PROFILE 


Inches 


Figure  164  TYPICAL  STATIC  PRESSURE  PROFILE 
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P  -  Static  Pressure  -  psia 


Figure  165 


TYPICAL  STATIC  PRESSURE  PROFILE  CORRECTION 
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Figure  166  THEORY/DATA  TOTAL  TEMPERATURE  PROFILE 

COMPARISON 
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Figure  167  THEORY/DATA  TOTAL  TEMPERATURE  PROFILE 

COMPARISON 


Figure  168  TYPICAL  TOTAL  TEMPERATURE  PROFILE  DATA 


TABLE  I 


MODEL  C  (RAMP  ONLY) 
PROFILE  ANALYSIS  SUMMARY 


RUN 

M0 

STA 

Ree 

6* 

0 

n 

Me 

(FT) _1 

(IN) 

(IN) 

22 

10.55 

1 

2.31X106 

.146 

.00560 

.960 

9.06 

10 

10.55 

2 

2.34X10® 

.173 

.00698 

3.08 

9.02 

14 

10.55 

3 

2.29X106 

.214 

.00941 

6.56 

9.00 

15 

10.55 

3 

2.62X106 

.214 

.00911 

6.22 

9.01 

20 

10.55 

1 

2.31X106 

.183 

.00770 

3.79 

9.07 

19 

10.55 

2 

2 .33X10® 

.262 

.01177 

10.0 

8.97 

21 

10.55 

2 

2.24X10® 

.242 

9.30 

9,00 

18** 

10.55 

3 

1.51X10® 

.322 

.01470 

11.2 

8.06 

23 

10.40 

1 

1.00X106 

.231 

.00870 

.942 

8.80 

11 

10.40 

2 

7.56X105 

.323 

.01270 

1.26 

8.42 

13 

10.40 

3 

8.37X105 

.324 

.01355 

1.24 

8.58 

71 

7.34 

1 

2.94X10® 

.0943 

.00691 

3.96 

6.50 

76 

7.34 

2 

3.31X10® 

.162 

.0126 

8.34 

6.60 

74 

7.34 

3 

2 .66X10® 

.220 

.0176 

7.95 

6.50 

69 

7.34 

1 

3.76X106 

.149 

.0106 

8.14 

6.54 

68 

7.34 

2 

3.22X10® 

.221 

.0166 

7.42 

6.58 

73 

7.34 

3 

3.60X10® 

.283 

.0212 

6.96 

6.55 

75 

7.34 

3 

3.29X10® 

.272 

.0210 

6.82 

6.55 

70*** 

7.30 

1 

4.94X105 

.207 

.0108 

.640 

6.20 

67 

7.30 

2 

4.97X105 

.201 

.0152 

1.14 

6.12 

72 

7.30 

3 

5.86X105 

.261 

.0194 

1.08L 

6.12 

2.08 


L  Lower  portion  of  boundary  layer 

U  Upper  portion  of  boundary  layer 

**  Survey  did  not  extend  to  boundary  layer  edge 


***  Survey  indicated  boundary  layer  separation 


TABLE  II 


MODEL  C 

PROFILE  ANALYSIS  SUMMARY 


RUN 

Mo 

STm 

ec 

Ree  , 

6* 

e 

n 

Me 

(DEGREES) 

(FT)-1 

(IN) 

(IN) 

48 

10.55 

3 

5 

4.71X106 

.143 

.01219 

10.4 

6.20 

49 

10.55 

3 

10 

4.96X106 

.077 

.01196 

11.3 

4.40 

50 

10.40 

3 

1.5 

1.03X106 

.218 

.01262 

3.7 

7.28 

51 

10.40 

3 

5 

1.70X10® 

.109 

.00848 

11.1 

6.34 

58 

7.34 

2 

5 

6 . 28X10° 

.109 

.01246 

5.84 

5.04 

61 

7.30 

2 

2 

9.14X105 

.140 

.01249 

1.24 

5.50 

62 

7.30 

2 

8 

1.12X10® 

.044 

.00729 

20.0** 

-- 

**  Outer  edge  of  boundary  layer 


TABLE  III 


MODEL  E-3 


PROFILE 

ANALYSIS 

SUMMARY 

RUN 

M0 

STA 

TRIP 

ec 

Ree 

8* 

e 

n 

Me 

(DEGREES) 

(FT)-1 

(IN) 

(IN) 

5 

10.40 

1 

-- 

6.58X105 

.219 

.00999 

1.13 

8.30 

6 

10.40 

1 

— 

6.08X105 

.210 

.0104 

1.18 

7.97 

8 

10.40 

2 

-- 

6.38X105 

.213 

.0116 

1.05 

7.30 

13 

10.40 

3 

— 

9.48X105 

.165 

.01105 

2.65 

7.22 

4 

10.55 

1 

#1 

-- 

2.09X106 

.141 

.00581 

1.27 

8.79 

10 

10.56 

2 

#1 

-- 

2.03X106 

.130 

.00707 

3.58 

7.90 

46 

10.56 

2 

#2 

-- 

1.87X106 

.146 

.00961 

3.45 

7.72 

11 

10.56 

3 

#1 

— 

2.65X106 

.151 

.0103 

11.5 

7.36 

45 

10.56 

3 

#2 

— 

2.33X10® 

.170 

.0125 

13.6 

7.27 

7 

10.56 

1 

— 

1.96X106 

.141 

.00607 

1.14 

8.50 

9 

10.55 

2 

-- 

2.31X106 

.133 

.00679 

3.50 

8.14 

12 

10.56 

3 

— 

2.56X106 

.156 

.0110 

11.8 

7.35 

29 

10.40 

3 

-1 

7.12X105 

.166 

.0121 

4.04 

6.68 

30 

10.55 

3 

U2 

-1 

2 . 18X10® 

.161 

.0118 

12.5 

6.96 

50 

10.56 

2 

* 2 

14.9 

1 .91X10® 

.141 

.0074 

9.17 

7.80 

38 

10.56 

3 

n 2 

15.1 

2.77X106 

.175 

.0115 

13.5 

7.57 

94 

7.30 

1 

— 

1.79X105 

.238 

.0205 

1.23 

5.56 

102 

7.29 

3 

— 

2.65X105 

.268 

.0264 

1.25 

5.41 

57 

7.30 

1 

— 

2.56X105 

.197 

.0176 

1.23 

5.22 

60 

7.30 

3 

— 

4 .21X10® 

.150 

.0157 

4.08 

5.36 

56 

7.34 

1 

— 

2.05X10® 

.108 

.0100 

1.09 

6.51 

54 

7.34 

2 

— 

2.35X10® 

.113 

.0111 

8.87 

6.27 

58 

7.34 

3 

— 

2.61X10® 

.123 

.0140 

11.3 

5.63 

100 

7.30 

2 

8.0 

2.21X105 

.239 

.0279 

1.42 

5.37 

79 

7.30 

3 

3.0 

2 . 64X10® 

.345 

.0308 

.907** 

4.53 

101 

7.30 

3 

8.2 

-- 

.498 

.0543 

Separated 

3.66 

78 

7.30 

3 

15.0 

3.17X105 

.180 

.0183 

1.29l 

4.71 

2.94u 

82 

7.34 

2 

15.0 

2.17X10® 

.112 

.0109 

9.40 

6.04 

92 

7.34 

3 

3.0 

2.36X10® 

.141 

.0163 

9.83 

5.63 

76 

7.34 

3 

15.0 

2.46X106 

.131 

.0143 

9.63 

5.65 

L  Lower  portion  of  boundary  layer 

U  Upper  portion  of  boundary  layer 

**  Survey  indicated  boundary  layer  separation 
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